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Observations  of  Vertically  Propagating,  Equatorial ly-trapped 
Waves  in  the  Deep  Western  Incian  Oc<?an 

ABSTRACT 

An  hypothesis  of  equatorial ly-trapped  waves  is  found  to  be 
consistent  with  time  series  of  vertical  profiles  of  horizontal  velocity 
and  CTD  data  from  the  western  Indian  Ocean.  The  profiles  were 
collected  using  an  acoustic  dropsonde,  the  White  Horse,  along  the  53°E 
meridian.  The  temporal  coverage  is  a  month-long  period  spanning  the 
onset  of  the  southwest  monsoon  in  1976;  the  latitudinal  coverage  is 
3/4°S  to  5°N. 

To  examine  the  composition  of  the  velocity  field,  a  WKB  stretching 
procedure  was  applied  to  the  depth  and  a  WK8  normalizing  procedure  to 
the  velocity.  An  average  buoyancy  frequency  profile  was  used  to  reduce 
each  velocity  profile  to  the  form  which  would  have  been  observed  in  a 
uniformly  stratified  ocean.  Removal  of  the  variations  of  amplitude  and 
scale  resulting  from  the  background  density  stratification  makes  it 
possible  to  discuss  vertical  wavenumbers  throughout  the  temporal  and 
spatial  array,  rather  than  merely  local  vertical  wavenumbers. 

Autospectral  estimates  reveal  equatorial  intensification  that 
varies  with  vertical  wavenumber.  The  hypothesis  is  that 
equatorial ly-trapped  waves  of  one  kind  or  another  are  the  basis  for  the 
observed  motion.  To  examine  vertical  propagation,  dropped  lagged 
coherences  have  been  computed.  Pairs  of  drops  separa.ed  by  the  same 


temporal  lag  were  grouped  together  to  estimate  cross-spectra.  The 
results  indicate  the  presence  of  a  mixed  Ross  by-gravity  wave  of  60-77 
day  period  with  phase  propagation  downward  at  a  vertical  wavelength  of 
1200  in  the  stretched  coordinate,  which  is  equivalent  to  approximately 
three  wavelengths  between  the  thermocline  and  the  bottom.  Vertical 
propagation  is  indicated  in  other  wavenumber  bands  as  veil,  also 
corresponding  to  autospectral  peaks. 

If  linear  equatorial  wave  theory  is  applicable,  zonal  wavelengths 
can  be  determined  from  the  dispersion  relation,  and  possible  forcing 
mechanisms  can  then  be  examined.  Linear  theory  gives  a  zonal 
wavelength  of  300-400  km  for  the  1200  stretched  meter  (sm) 
oscillation.  The  rms  velocity  is  the  same  order  as  the  phase  speed, 
oowever,  so  that  non  linearities  may  not  be  completely  negligible. 
Indications  are  of  multiple  processes  within  the  same  wavenumber  band, 
a  longer-period  Kelvin  or  Rossby  wave  coexisting  with  the 
above-mentioned  mixed  Ross  by-gravity  wave. 
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Chapter  I.  Introduction  and  Theory 

Section  A.  Introduction 

A  unicue  set  of  data  was  collected  during  May  and  June  of  1976: 
vertical  profiles  of  horizontal  velocity  between  the  surface  and  the 
bottom  of  the  western  equatorial  Indian  Ocean.  The  discovery  that  the 
equatorial  current  structure  was  energetic  to  great  depths  was  the 
initial  rtsult  from  these  observations.  Instead  of  merely  adding  to 
our  knowledge  of  the  generation  and  temporal  structure  of  the 
undercurrent  in  the  Indian  Ocean,  Luyten  and  Swallow  (1976)  documented 
the  first  evidence  of  deep  vertical  structure  in  both  horizontal 
velocity  components.  Previous  observations  at  the  equator  had 
indicated  that  the  zonal  velocity  field  had  variability  below  the 
undercurrent,  but  profile  measurements  had  only  penetrated  to  a  maximum 
depth  of  1500m  (Rual,  1969).  The  type  of  measurement,  i.e.,  direct  and 
to  the  bottom,  possible  with  the  White  Horse,  a  freely  falling  acoustic 
dropsonde  developed  at  Woods  Hole,  represented  a  technological  and 
logistical  breakthrough. 

A. 1.  Brief  description  of  the  data.  The  White  Horse  was  deployed 
forty-one  times  at  seven  stations  in  the  region  between  50’30'E  and 
5 3 "E ,  3/4*S  and  5*N.  The  data  collected  were  unique  in  that  they 
provided  the  first  detailed  snapshots  of  the  equatorial  water  column 
and  because  of  the  complex  vertical  structure  they  revealed.  They 
remain  unique  because,  although  the  instrument  and  others  resembling  it 
have  been  ised  since  then,  there  is  still  no  comparable  time  series,  in 


quality  and  in  latitudinal  coverage  over  time. 

Figure  1-1  shows  the  location  of  the  profiling  stations.  The 
three-dimensional  position  of  the  instrument  was  determined 
approximately  every  twenty  meters  of  depth,  in  a  manner  described 
below,  and  horizontal  velocities  were  calculated  from  the  positions. 
The  resulting  vertical  profiles  of  horizontal  velocity  between  the 
surface  and  the  bottom  are  non-uniformly  separated  in  time.  (See 
Figure  1-2,  ship's  track.) 

Figure  1-3  displays  most  of  the  velocity  data.  The  drops  were 
primarily  along  the  53°E  meridian,  and  ranged  latitudinally  from  3/4*S 
to  5"N.  The  temporal  coverage  is  from  May  14  to  June  18,  1976.  The 
longest  time  series  at  a  single  location  is  at  the  equator,  with  ten 
profiles  at  0*,  53°E.  The  figure  shows  perspective  plots  of  the 
velocity  vectors  from  the  surface  to  the  bottom,  and  reveals  that  the 
vertical  structure  is  complex  at  the  equator  and  becomes  simpler  away 
from  it.  The  velocity  vectors  rotate  with  depth,  but  not  with  any 
regularity. 

Part  of  the  instrumentation  package  of  the  White  Horse  is  a  Neil 
Brown  conductivity,  temperature  and  depth  recorder  (CTO),  so  concurrent 
density  profiles  were  taken  from  the  surface  to  the  bottom.  The 
profiler  data  are  supplemented  with  data  from  current  meters  and 
temperature /pressure  recorders,  moored  in  the  same  region  where  the 
White  Horse  was  deployed.  Additional  data  were  collected  during  two 
INDEX  cruises  in  the  spring  and  summer  of  1979.  Reference  to  results 
from  those  experiments  will  be  made  where  appropriate. 
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FIGURE  1-1 

8athymetr1c  map  of  western  equatorial  Indian  Ocean, 
showing  location  of  1976  White  Horse  profiling  stations. 
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Ordinate  Indicates  latitude  (and  the  longitudinal 
exception,  net  5  at  50*30' E),  and  White  Horse  profile 
stations  are  labelled  with  the  net  number  followed  by 
two  sequential  digits. 
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Several  features  are  obvious  in  the  data.  The  vertical  structure  of 
the  velocity  field  at  the  equator  is  complex,  whereas  the  vertical 
structure  at  5*N  is  similar  to  a  first  baroclinic  mode,  i.e.,  the  complex 
vertical  structure  is  equatorial ly  confined.  The  large  scale  structures 
in  the  zonal  velocity  profiles  persist  with  time,  indicating  that 
long-period  processes  dominate  the  east-west  component.  The  velocity 
components,  u  and  v,  are  of  the  same  order,  rather  than  being 
overwhelming ly  zonal,  as  was  found  in  the  western  Pacific,  for  example 
(Eriksen,  1981).  A  dominance  of  clockwise  or  counterclockwise  rotation 
occurs  in  bands,  and  is  not  uniform  from  the  surface  to  the  bottom. 

When  the  velocity  profiles  are  examined  in  componen  .  form,  rather 
than  in  the  vector  form  shown  in  Figure  1-3,  a  dichotomy  can  be  seen  in 
both  the  vertical  and  the  temporal  scales  of  the  motion.  The  meridional 
component  displays  shorter  scales  in  both  dimensions.  Figure  I-4a  shows 
the  east-west  velocity  component  at  the  equator  in  original  form.  Time 
is  increasing  to  the  right.  The  first  impression  is  one  of  a  wavy 
structure,  with  higher  amplitudes  and  shorter  vertical  scales  in  the 
upper  part  of  the  water  column.  The  same  pattern  holds  for  the 
meridional  component  (Figure  I— 4b ) .  A  background  dependence  upon  the 
buoyancy  frequency,  which  is  indicated  very  clearly  in  the  current  meter 
analysis  done  by  tuyten  (1982),  seems  to  be  affecting  both  velocity 
amplitudes  and  vertical  length  scales. 

In  an  attempt  to  remove  this  background  dependence  and  to  reduce  each 
velocity  profile  to  the  form  which  might  have  been  observed  in  a 
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FIGURE  1-4 

Zonal  (a)  and  meridional  (b)  velocity  profiles  at  the 
equator,  53  E.  Time  Increases  to  the  right,  but  the 
spacing  between  profiles  Is  not  adjusted  for  the  actual 
time  Interval  between  them. 
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uniformly  stratified  ocean,  the  average  buoyancy  frequency  profile  was 
used  to  scale  both  velocity  and  depth.  The  "after"  picture  in  Figure 
1-5  shows  the  u-  and  v-components  after  this  WKB  stretching  and 
normalizing  procedure.  These  are  relabelled  u*  and  v*,  and  plotted 
against  z*.  The  oscillations  are  much  more  uniform,  except  in  the  top 
1400  scaled  meters,  which  correspond  to  the  top  275m  of  the  ocean.  The 
upper  portion  of  the  water  column  will  in  general  not  be  considered  in 
the  following  analysis. 

One  can  now  discuss  vertical  wavenumbers  throughout  the  temporal 
and  spatial  array  rather  than  merely  local  vertical  wavenumbers.  The 
data  in  this  form  are  the  inputs  to  the  spectral  analysis  discussed  in 
Chapter  1 1— C  and  [)  and  interpreted  in  Chapters  III  and  IV.  Variations 
of  the  cross-spectra  computed  from  the  profiles  at  a  particular 
vertical  wavenumber  should  then  depend  only  on  the  time  differences 
between  profiles.  An  analysis  can  be  done  using  all  lowerings 
separated  by  the  same  time  lag.  [Hayes  (1975)  analyzed  temperature 
profiles  in  this  fashion.]  For  each  temporal  lag  consicered,  coherence 
and  phase  can  be  examined  as  a  function  of  vertical  wavenumber,  m. 

One  can  expect  coherence  to  drop  as  m  increases,  since  the  data 
collection  method  operates  as  a  low-pass  filter  by  averaging  the 
velocity  vector  over  20-second  intervals  (Hendricks  and  Rodenbusch, 
1981),  and  since  smaller  wavelengths  should  be  less  correlated  as  the 
band-width  increases  (Munk  and  Phillips,  1968).  The  coincidence  of 
zero  phase  with  strong  coherence  would  indicate  no  phase  propagation  at 
that  wavenumber.  Phase  change  for  a  particular  wavenumber  with  longer 
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U*  (scm/ttc) 


V*  (scm/s»c) 


AFTER 

FIGURE  1-5 

The  same  zonal  (a)  and  meridional  (b)  velocity  profiles 
as  In  Figure  1-4,  after  a  WKB  stretching  and  normalizing 
procedure  was  applied. 
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temporal  lags,  from  which  vertical  propagation  can  be  inferred,  is 
discussed  in  Chapter  II-D.2. 

The  purpose  of  this  dissertation  is  to  examine  observations  of 
long-period  motion  in  the  equatorial  zone  and  to  infer  from  these 
observations  the  governing  dynamics.  The  primary  emphasis  will  be  on 
analysis  of  the  profiling  data.  The  hypothesis  is  that  long-period 
equatori ally-trapped  waves  of  one  kind  or  another  are  the  basis  for  the 
observed  motion.  There  is  a  presumption  against  the  existence  of 
vertical  modes,  based  on  the  arguments  of  Wunsch  (1977),  summarized  as 
follows.  Theoretically,  a  wave  of  long  period  propagates  energy  from 
the  surface  at  only  a  shallow  angle  to  the  horizontal  and  would  take 
several  years  to  reach  the  ocean  floor.  Dissipation  would 
substantially  reduce  the  energy  of  the  wave  before  normal  modes  could 
be  set  up,  even  if  the  wave  were  perfectly  reflected  frcm  the  bottom 
and  from  meridional  boundaries.  On  the  other  hand,  vertical 
propagation  is  not  visually  evident.  If  there  were  a  dominant 
semi-annual  signal,  for  example,  a  60°  phase  change  could  be  observed 
over  the  month-long  time  domain  involved,  yet  this  is  not  the  case. 

It  is  possible  to  imagine  generation  mechanisms  for  the  White  Horse 
profiles  other  than  vertically-propagating  linear  equatorial  waves,  but 
one  feature  of  the  data  that  must  be  explained  is  equatorial 
intensification,  a  highly  characteristic  feature  of  such  waves.  Any 
discrepancies  between  linear  theory  and  observations  which  have  been 
clearly  and  accurately  identified  are  pursued  as  far  as  possible  here 
and  will  motivate  further  theoretical  research. 
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A. 2.  Organization.  The  basic  approach  outlined  above  leads  to  the 
following  structure  for  the  dissertation.  The  remincer  of  Chapter  I 
presents  a  review  of  linear  equatorial  wave  theory.  Chapter  II 
describes  the  methods  used  in  analyzing  the  data,  and  presents  some 
interpretative  guidelines.  Because  of  their  differences  in  dominant 
frequencies  and  vertical  wavenumbers,  the  velocity  components  are 
handled  in  separate  chapters.  Chapter  III  presents  the  results  from 
analysis  of  the  zonal  velocity  component.  Chapter  IV  the  meridional. 
Chapter  V  contains  a  comparison  with  other  results  and  recapitulates 
the  primary  conclusions  from  this  research. 
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Section  8.  Linear  Equatorial  Wave  Theory  -  Review 

B.l.  General  linear  theory  and  solutions.  The  linear  Boussinesq 
equations  for  inviscid,  incompressible,  unforced  rotating  flow  on  an 
unbounded  equatorial  beta-plane  are: 

ut  -  eyv  =  -px  (la) 

vt  +  syu  =  -py  (lb) 


o  =  -p_  -  asL 


Z  P, 


(1 C) 


pt  +  WP2  =  0 


UX  +  Vy  +  WZ 


(ld) 

(le) 


where  the  subscript  denotes  differentiation,  p  is  reduced  perturbation 

pressure  (p  =  p'pQ),  and  p  =  p  +  p ' ;  pQ  is  the  constant  reference 

density  and  o  is  the  observed  density  field,  composed  of  mean  (p) 

and  fluctuating  (p1)  components.  By  cross-differentiation  and 

using  R2  =  -  ^  ,  equations  lc  and  Id  reduce  to  p  .  +  R2w  =  0. 

po  32  zt 

Fo'lowing  the  notation  of  Eriksen  (1982),  let 


u 

— 

U  (y/y*)G'(z) 
nu,,r  mv  ' 

V 

,vn(y/ys)G;(z) 

p 

A  i (kx-wt) 
rtnme 

w 

.  1. 

(2) 


where  cm  =  Jgbm  [=  NQ/ |m )  ]  is  a  separation  constant,  with  units  of 


Vi 


speed;  hm  is  defined  as  the  equivalent  depth;  NQ  is  a  reference  buoyancy 
frequency;  y*m  =  [c  /b]*^  is  a  characteristic  meridional  length  scale, 
the  equatorial  radius  of  deformation;  and  5  is  the  vertical  displace¬ 
ment  of  isopycnals,  measured  upward  (see  Chapter  1 1— B  for  derivation 
from  the  conservation  of  mass  equation,  Equation-Id).  The  dependent 
variables  are  assumed  to  be  simply  periodic  in  x  and  t,  and  some 
as-yet-undef ined  functions  of  y  [Un,Vn  or  Pp(y/y*)]  and  z  [G  or 
3'(z)]. 


B.l.a.  Vertical  structure  equation.  The  equations  are  separable, 
ind  reduce  to  a  vertical  and  a  meridional  structure  equation.  The 
horizontal  structure  equation  is  similar  to  that  derived  in  shallow 
water  wave  theory  from  the  Laplace  Tidal  Equations,  where  here  the 
depth  parameter  is  hffl  (*  equivalent  depth)  rather  than  H  (=  total 
water  depth).  The  vertical  structure  equation  is 


+  E(z)6  _ 


zz 


=  o  , 


(3) 


m 


which  admits  first-order  WKB  solutions  of  the  form 

GU)  *  -jf-i/2  /  r;  d2>  ■ 

[?-]  J  m 


where  N(z)/cm  =  (ml  is  the  local  vertical  wavenumber  (Philander,  1978) 
and  the  choice  of  sign  depends  upon  the  direction  of  propagation  [plus 
(minus)  indicates  phase  propagation  upward  (downward)].  If  we  let 
m*  =  NQ/cm  then  the  approximate  solution  becomes 


G(z)  * 


!m*|1/2  (  ff(z) 


|  — —  |  exp(sgn  m*i  j  |m*|  dz)  . 


(4) 
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A  new 

variable  z* 

can  then  be  defined: 

dz* 

=  N(z)/Nq  dz. 

Thus 

z* 

=  / N(z)/N0  dz 

and 

G(z) 

=  [yNUH^GU*); 

G*  (z 

)  =  sgn  m  [N(z)/N0]1/2G‘ (z*). 

Thus 

u,  V,  p  OC 

[N{z)/N0]1/2G'(z*);  wj  <£  [N0/N(z)]1/2G(z*) 

Th«‘  WKB-scaled  variables  are: 


u*l  f  N 

o 

»'★  _  _ 

'  ‘  N(z) 


WKB 

Scaling 


“  -1 

"  -  Mil 

[  J*j  '  No 


i  L  n°  j  m  ’  _ 

with  a  scaled  vertical  wavenumber  |m*(z*)|»  NQ/  J ghm,  which  is  no 
longer  dependent  upon  the  local  buoyancy  frequency,  but  only  upon  the 


(8) 


15 


where  n  -  0,1,2,....  A  "turning  latitude",  yT,  can  then  be  defined, 
where  the  form  of  the  equation  changes  from  elliptic,  with  oscillatory 
solutions,  to  hyperbolic,  with  solutions  which  decay  exponential ly 
uoleward. 


yT 


2 


(2n  +  1)  . 


(9) 


This  critical  latitude  will  be  discussed  qual itati vely  later,  with 
respect  to  particular  meridional  modes. 

With  the  transformation  of  variables  above,  the  dispersion  relation 
(Equation-I-8)  can  be  expressed  as  a  function  of  vertical  wavenumber: 


5!^!  _  k2  -  lii  =  4^  (2n  +  1)  ,  (10) 

and  the  turning  latitude  is  inversely  proportional  to  the  square  root 
of  the  vertical  wavenumber: 


yT  = 


l  |">l  8 


(2n+l)N0  1/2 

- 


(ID 


The  meridional  structure  equation  reduces  to  canonical  form  with 
any  of  the  equivalent  substitutions: 


n 


[7F] 


r3  1ml  ! 

[_^r3 


y/y*  , 


N0  1/2 

where  y*  *  [— j-|  ]  is  the  characteristic  meridional  length  scale, 

also  known  as  the  equatorial  Rossby  radius  of  deformation.  Its 
solutions  are  Hermite  functions: 
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vn(y/y*)  =  An /n(y/y*) 

exp(-y2/2y*2)Hn(y/y*) 


=  A. 


n  _  "(  2^  nl  *  w  1/2)  1/2  * 

where  Hn(y/y*)  are  the  Hermite  polynomials  of  order  n.  [HQ=1,  H1=2n, 
H2=4n2-2,  H2=8n^-12n,  ....]  The  Hermite  functions  form  a  complete 
orthogonal  set,  and  with  the  above  normal ization  factor. 


(12) 


OQ 

/ 


y m /'n  dn  =  smn  ’  the  Kronecker  delta. 


The  corresponding  solutions  for  Un  and  Pn,  after  substitution 
into  the  original  set  of  equations  (1-1),  are: 


u n(y/y*) 


Pn(y/y*)  *  a 


r  n  ^2 

A  r  <s»  t«-i 

n  |  o  +  s 

r  „  1/2 

[ifLJVi 

n  I  o  +  s 


,n+l 


1/2 


<¥>  fn*l 


o  -  S 
1/2 


(!!r)  KU 

a  -  1 


(13a) 


(13b) 


where,  for  simplicity  of  presentation,  the  meridional  structure  equation 

J 1/2 

has  been  non-dimensional ized  using  a  =  — 5:1 — m  and  s 

Ki177 


i  uv* uurc  ci 

-m 


k. 


as  well  as  n  =  y/y* 


“  • 
6  I  m| 

=  iVj 


1/2 


y,  to  reduce  it  to  a  canonical  form: 


B.2.  Particular  solutions  and  discussion.  The  solutions  to  this 
equation  are  of  three  types,  each  corresponding  to  a  root  of  the  cubic 
equation  for  frequency,  the  dispersion  relation  (1-10).  These  solutions 
will  be  discussed  following  development  of  a  fourth  type  of  solution  to 


r 
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equation  (1-1). 

B.2.a.  Kelvin  waves.  An  equator i ally-trapped  wave  type  which  does 
not  arise  directly  out  of  the  above  analysis  (since  V  =  0),  is  the 
Kelvin  wave.  By  convention,  it  is  ascribed  a  meridional  mode  number  of 
-1,  since  for  n  *  -1, 
ir  to  to  Imj 

particular  wave  solutions  for  the  Kelvin  wave  are 


is  a  solution  to  the  dispersion  relation.  The 


a-y2/2y*2 


u_1(y/y*)  =  P.^y/y*)  =  A_x  -c  r/T 


so  that 


and 


-lr 


u_1(x,y/y*,z,t)  =  cm~  P_1(x,y/y*,z,t) 


-1, 


^(x.y/y*. z.t)  =  cffl  P_^(x,y/y*,t)  6 (z) 


=  P_1(x,y/y*,z,t) 


.-1. 


=  -i  rrf  u_1(x,y/y*,z,t)  . 

The  meridional  structures  of  the  u  and  p  fields  are  Gaussian  in  latitude, 
identical  except  for  an  amplitude  factor.  The  displacement  field  is  in 
quadrature  with  both  of  them,  leading  u  by  n/2  with  depth.  The  Kelvin 
wave  zonal  phase  velocity  (w/k)  is  thus  cm  =  =  N0/|m|  ,  directed 

toward  the  east,  and  exactly  equal  to  the  zonal  group  velocity  (Wak). 
The  vertical  phase  velocity  (u»/m)  depends  directly  upon  the  sign  of  the 
vertical  wavenumber  (upward  if  m  >  0).  Vertical  group  velocity  is 
equal  in  magnitude  and  opposite  in  sign  to  the  phase  velocity. 
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B.2.b.  Mixed  Rossby-gravity  waves.  When  n  =  0,  the  dispersion 
relation  ( I — 10 )  can  be  factored.  The  solution  for  which  th?  associated 
zonal  velocity  field  is  not  singular  is 

a*) 

The  wave  behaves  as  a  gravity  wave  at  high  frequencies  and  as  a  Rossby 
wave  at  low  frequencies,  and  consequently  is  known  as  a  mixed 
Rossby-gravity  wave  (or,  by  meteorologists,  sometimes  called  a  Yanai 
wave).  The  particular  wave  solutions  are 


Vo  =  flofo=T74ex>,(-*2,2y*2) 

IT 

A„ 


U  =  P  =  - 
0  0 


/r 


0  _ 1 

O-S 


>  -  exp (-//2y*2)  . 

ir  0 

The  zonal /meridional  kinetic  energy  ratio  at  a  given  latitude  is 
2  2  2  2  2 

thus  u0  'Vq  *  —  •  The  zonal  phase  velocity  is  given  by 

No 

w2N0/(u>2|ml  -  BNQ),  directed  toward  the  west  at  low  frequencies,  toward 
the  east  at  high  frequencies.  The  zonal  group  velocity  is  slower,  and 
always  directed  to  the  east:  w2N0/(w2|ml  +  bNq) .  The  direction  of 
the  vertical  phase  velocity  is  dependent  upon  the  sign  of  the  vertical 
wavenumber:  (sgn  m)  u>3/(kN0w  +  eNQ).  The  vertical  group  velocity  is 
slower  and  in  the  opposite  direction:  (-sgn  m)u>3/U2/mf  +BN0). 

The  ratio  of  zonal  to  vertical  group  velocity  is  easily  seen  to  be 
(-sgn  m)  N 0/«,  for  both  the  Kelvin  and  the  mixed  waves.  Since  the 


frequencies  of  particular  interest  in  the  equatorial  oceans  are  always 
much  less  than  NQ  (~  one  cycle  per  hour),  the  ratio  indicates  that 
the  path  along  which  energy  propagates  for  these  equatorial ly-trapped 
waves  will  always  be  at  only  a  shallow  angle  to  the  horizontal , 
shallower  for  longer  periods. 

The  characteristic  meridional  scale,  or  LR,  the  equatorial 
Rossby  radius  of  deformation,  defines  the  e-folding  distance  from  the 
equator  for  the  kinetic  energy  in  the  zonal  velocity  component  of  a 
Kelvin  wave  and  for  the  meridional  component  of  a  mixed  Rossby-gravity 
wave.  Note,  from  Table  I— 1,  that  LR  is  less  than  the  distance 
between  the  equator  and  the  nearest  off-equatorial  profiling  station 
(3/4°N  or  S)  for  vertical  wavelengths  shorter  than  514  stretched 
meters,  so  that  the  meridional  structure  at  such  wavelengths  will  be 
more  difficult  to  observe  with  the  data  analyzed  here. 

B.2.C.  Rossby  and  inertia-gravity  waves.  There  are  two  classes  of 
equatorial ly-trapped  waves  for  meridional  mode  numbers  equal  to  1  or 
larger,  separated  by  a  frequency  range  which  can  be  filled  only  by  the 
Kelvin  and  mixed  waves  discussed  above.  At  low  frequencies  (a  <  1), 
these  are  the  planetary  or  Rossby  waves;  at  high  frequencies  (o  >  1), 
these  are  the  inertia-gravity  waves  (similar  to  mid-latitude  internal 
gravity  waves  but  trapped  by  the  beta-effect  to  the  equator).  The 
minimum  period  at  which  Rossby  waves  of  specific  vertical  wavelengths 
can  exist,  and  the  corresponding  maximum  for  inertia-gravity  waves,  can 
be  determined  from  the  dispersion  relation  (1-10),  by  solving  for  the 
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TABLE  1-1 


2w 

N 

,  0 

X  . 

z* 

"^7 

hm=  X 

z 

m  g 

(sm) 

(sm)"1 

l 

(cm) 

3600 

1.7  x 

10-3 

12.18 

1800 

3.04 

1200 

5.2  x 

10-3 

1.35 

900 

7.0 

.76 

720 

8.7 

.49 

600 

1.0  x 

10-2 

.34 

514 

.25 

450 

1.4 

.19 

400 

1.6 

.15 

360 

.12 

327 

.10 

300 

2.1 

.08 

277 

257 

240 

2.6 

.05 

225 

212 

200 

3.1  x 

10-2 

.04 

Length  Scale 

Time  Scale 

v 

[s=*2h] 

lr 

[c=2ir] 

T 

(cm/sec) 

(km) 

(days) 

109 

218 

2.3 

55 

154 

3.3 

36 

126 

4.0 

27 

109 

4.6 

22 

98 

5.2 

18 

89 

5.7 

16 

83 

6.1 

14 

77 

6.5 

12 

73 

6.9 

11 

69 

7.3 

10 

66 

7.7 

9 

63 

8.0 

61 

8.3 

58 

8.7 

7 

56 

9.0 

55 

9.3 

53 

9.5 

6 

51 

9.8 

,  ll/2 

Jj 

(13.24579  x  xz 

[m])1^2  km 

m2g 


(.93988  xx22  [km])  cm 


N0  -  1.091464  cph 
*  26. 195  cpd 

-  164.6  d-1 

-  1.9049  10'3  s-1 


9g»  =  978  cm  sec-2 

80»  =  2.289  10"°  cm"1  s"1 
-  1.9777  10-6  m-l  d-l 
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zonal  wavenumber  and  examining  the  term  enclosed  under  the  radical: 

~  2  2  2l 

k  =  -  £-  ±  -  ilsL  (2n+l)  +  2-S-  .  (16) 

4*>  NnZ 

v  0  o-* 

For  n=l,  the  two  values  of  u  which  make  the  expression  in  brackets 
equal  to  zero  can  be  determined: 


The  upper  factor  defines  the  maximum  inertia-gravity  wave  period  (for 
the  first  meridional  mode);  the  lower  defines  the  minimum  period  for 
the  first  mode  Rossby  wave. 


These  maximum  and  minimum  periods  are  shown  in  Figure  1-6  (connected 
by  a  dotted  line),  where  the  dispersion  relation  for  xz  =  1200sm  is 
plotted  as  a  function  of  frequency  and  zonal  wavelength.  The  dotted 
line  is  defined  by: 


which  is  where  the  zonal  group  velocity  is  exactly  zero  (since  the  term 
under  the  radical  in  Equation  1-16  is  identically  zero).  The  curves 


792  km 


i  i - 1 — 

-464  -792  -2700 


Figure  1-6 

A  cut  of  the  three-dimensional  dispersion  relation 
(Equation  1-10)  at  the  vertical  wavenumber  corresponding 
to  a  wavelength  of  1200sm.  The  curves  are  labelled  by 
wave  type,  and  critical  periods  and  zonal  wavelengths  are 
indicated. 
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corresponding  to  Kelvin  and  mixed  Ross by-gravity  waves  are  labelled,  as 
well  as  those  for  the  Rossby  and  inertia-gravity  waves.  The  phase 
velocity  of  Rossby  waves  is  seen  to  be  always  westward,  and  the  zonal 
group  velocity  is  westward  for  long  zonal  wavelengths  and  eastward  for 
short  wavelengths.  The  dotted  line  distinguishes  the  inertia-gravity 
waves  on  the  right,  which  propagate  energy  to  the  east,  and  those  on 
the  left  which  propagate  energy  to  the  west. 

The  general  forms  for  zonal  and  vertical  group  velocity  are: 
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x 
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2ukN  2  +  bN  2 

o _ o _ 

3<u2m2  -  k2NQ2  -  |m| bN( 2n+l) 


and 


-  S  -  (s9n  m) 


am 


«uBN02(2n+l)  -  2|m|(u3 
3«V  -  k2N  2  -  |m(BN(2n+l) 


Since  the  denominators  are  the  same,  the  ratio  of  the  zonal  group 
velocity  to  the  vertical  group  velocity  is  simply 
x  o  2  +  eN  2 


2u.kNn 

-3-  =  (sgn  m)  - * 

c„z  <uBN0(2n+l)  -  2lm|uJ 


c  x  (sgn  m)N 

v/hich  reduces  to  =  -  -^^+Yj "  the  ^ow  frequency  limit  for  Rossby 
c 

waves.  3 

The  Gaussian  term  in  the  general  wave  solutions  (Equations  1-12, 

13a,  13b),  which  dominates  the  meridional  structure  of  the  zonal 
velocity  of  the  Kelvin  waves  and  the  meridional  velocity  of  the  mixed 
waves,  acts  as  a  modifying  factor  rather  than  an  envelope  for  the 


Rossby  and  inertia-gravity  modes.  Their  structure  is  better  described 
in  terms  of  the  turning  latitude  (Equation  I— 11 )  than  in  terms  of  the 
e-iolding  distance.  The  meridional  mode  number  defines  the  number  of 
nodes  in  the  V  field  between  the  equator  and  the  poles,  and  thus 
functions  in  the  place  of  a  meridional  wavenumber.  The  turning 
la.itude  defines  the  location  where  the  wave  solutions  begin  to  decay 
poleward. 

For  inertia-gravity  waves  of  a  particular  meridional  mode  number, 
the  longer  the  period  the  more  tightly  trapped  to  the  equator  the  waves 
become,  since  w  =  f  defines  the  turning  latitude.  Equation  1-8  can 
be  solved  for  (ghm)l/2. 


(9hJ 


-u(2n+l)  ±  [fs2(2n+l)2  +  4(—  +  k2)w2]1/2 


2(k  ♦  ) 

ID 


and  the  inertia-gravity  solution  corresponds  to  the  choice  of  the  plus 

2 

sign.  This  equation  reduces  to  u>  /e(2n+l)  for  large  values  of  n,  so 
that  Equation  1-9  becomes 


On  the  other  hand,  the  number  of  Rossby  modes  which  can  be  excited 
increases  as  the  frequency  decreases.  If  the  equivalent  depth  (h") 
is  small  (i.e.,  if  the  vertical  length  scale  is  short),  the  waves  are 
equatorial ly-trapped,  since  the  turning  latitude  becomes  further  from 
thr  equator  as  the  equivalent  depth  increases.  For  large  values  of 
hm  or,  equivalently,  long  vertical  length  scales,  the  Rossby  waves 
cannot  be  adequately  described  on  an  equatorial  beta-plane.  These 
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limits  were  pointed  out  by  Philander  (1978). 

B.2.d.  Symmetry  and  antisymmetry.  The  meridional  symmetry  of  the 
problem  defined  by  the  different i a  1  equation  and  the  boundary 
conditions  implies  that  the  solution  will  be  composed  of  symmetric  or 
antisymmetric  functions.  For  the  meridional  velocity  field,  the 
meridional  mode  number  is  the  Hermite  polynomial  order  number.  It  has 
odd  or  even  symmetry  about  the  equator  depending  upon  the  mode  number, 
because  the  Hermite  polynomials  of  odd  order  are  sums  of  odd  powers  of 
latitude.  The  zonal  velocity  and  pressure  wave  fields  have  opposite 
symmetry  to  the  v  field,  since  they  are  composed  from  the  Hermite 
polynomials  of  adjacent  orders.  For  meridional  modes  of  odd  order, 
the  v  field  has  a  node  (zero)  at  the  equator  (i.e.,  is  antisymmetric 
about  the  equator)  and  the  u  and  p  fields  have  a  local  maximum  or 
minimum  at  the  equator  (i.e.,  are  symmetric  about  the  equ  itor  • .  The 
converse  holds  a  >  well. 


CHAPTER  II  -  Data  Acquisition  and  Reduction  Procedures 


The  methods  used  to  reduce  the  White  Horse  data  for  subsequent 
analysis  are  discussed  here  after  a  brief  description  of  the  instrument 
and  its  deployment.  The  data  recorded  by  the  instrument  are  first 
converted  to  velocity  components  and  the  horizontal  velocity  profiles 
are  then  scaled  using  the  buoyancy  profile  calculated  from  the  density 
data.  Vertical  displacement  is  also  calculated  from  the  density  data, 
and  similarly  scaled,  in  order  to  provide  additional  comparisons  to 
linear  wave  theory.  The  WKB  scaling  transformations  are  chosen  to 
provide  constant  coefficients  to  the  lowest  order  in  the  vertical 
structure  equation,  allowing  sinusoidal  solutions,  for  use  in  possible 
interpretation  of  vertical  wavenumber  analysis  of  the  data. 

Section  A.  Initial  Data  Processing 

A. 1.  Description  of  the  instrument  and  its  deployment.  The  White 
Horse  System  consists  of  an  acoustic  dropsonde,  approximately  two 
meters  long  and  containing  a  standard  Neil  Brown  CTD,  acoustic 
navigation  and  recording  equipment,  and  a  network  of  (initially)  three 
bottom-tethered  transponders.  Figure  1 1— 1  is  a  photograph  of  the 
instrument  being  lowered  over  the  side  of  a  research  vessel,  after 
which  it  is  released  and  falls  freely  to  the  bottom.  The  transponders 
have  previously  been  deployed  in  an  approximately  equilateral  triangle 
(on  the  order  of  4  km  per  side)  on  the  ocean  bottom.  As  the  dropsonde 
descends  through  the  water  column,  it  transmits  an  11.0-Khz  pulsed 
signal  every  20  seconds,  which  stimulates  a  return  signal  from  each  of 
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FIGURE  II-l 


Photograph  of  the  White  Horse  on  a  crane, 
prior  to  launching  at  a  profile  station. 
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the  bottom  transponders.  The  time  intervals  between  interrogation  and 
reply  are  stored  on  a  magnetic  tape  cassette  within  the  White  Horse. 

When  the  instrument  reaches  the  ocean  floor,  its  anchor  is  released 
by  one  of  three  mechanisms:  a  pre-set  clock,  an  acoustic  signal  or  a 
corrosible  magnesium  link  in  the  connecting  chain.  The  instrument  then 
rises  to  the  surface,  acoustically  tracking  its  trajectory  as  in 
descent.  An  acoustic  signal  to  the  AMF  release  is  the  normal  method 
and  reliance  on  the  corrosible  link,  the  second-order  backup,  has  never 
been  required  for  retrieval. 

An  acoustic  survey  is  performed  from  the  surface  to  determine  the 
depths  and  relative  positions  of  the  three  bottom  transponders.  The 
relative  positions  of  the  transponders  are  resolved  using  a  least 
squares  minimization  procedure,  from  an  overdetermined  (more  than  six 
interrogation  sites)  set  of  ranges  from  the  surface  to  the  transponders 
(vanderKulk,  1961).  Additional  information  collected  during  the  survey 
from  satellite  fixes  and  dead-reckoning  navigation  allows  the 
subsequent  determination  of  the  absolute  positions  of  the  transponders, 
and  thus  the  absolute  orientation  of  the  network,  to  within  ±5°.  Table 
1 1— 1  contains  survey  results  for  the  1976  White  Horse  Indian  Ocean 
networks. 

A. 2.  Conversion  of  travel  times  to  velocity  components.  The 
round-trip  travel  times  are  converted  to  three-dimensional  ranges  using 
the  local  sound  velocity  profile,  computed  from  the  CTD  data.  An 
iterative  triangulation  calculation  (described  in  Luyten,  Needell  and 
Thomson,  1982)  uses  the  ranges  and  transponder  positions  to  find  the 
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TABLE  II-l 

WHITE  HORSE  SURVEY  RESULTS 


Prof i  le 

Location 

Date  of 

Net  Coordinates 

Orientation 

Series 

(Beacon  1 

Survey 

beac  1:  x,y,z 

of  Net 

approx. ) 

1976 

beac  2:  x,y,z 
beac  3:  x,y,z 

*E  of  N 

500 

0°  N 

May  13 

0,  0,5049 

9 

50°30'E 

5291,  0,5046 

1551,3950,5052 

600 

0‘  N 

May  15 

0,  0,5049 

0 

53’  E 

4369,  0,5043 

1323,3228,5049 

0*  N 

June  15 

0,  0,5049 

0 

53°  E 

3620,  0,5043 

700 

0*45'N 

May  29 

0,  0,5054 

0 

53*  E 

2137,  0,5053 

1154,1715,5057 

800 

1  *  30 1 N 

May  30 

0,  0,5091 

-5 

53*  E 

2192,  0,5085 

1129,1834,5091 

900 

3*  N 

May  31 

0,  0,5099 

7 

53’  E 

3547,  0,5101 

1745,2161,5101 

1000 

5*  N 

June  1 

0,  0,5083 

5 

53’  E 

3792,  0,5083 

2457,2181,5085 

1100 

0’45'S 

June  5 

0,  0,4994 

0 

53*  E 

3880,  0,4991 

2667,3439,4993 

Net 

Dimensions 

(km) 


-31/2- 
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(x,y,z)  location  of  the  instrument  within  the  water  column.  The 
instrument’s  position  is  determined  every  20  seconds,  or  approximately 
every  20m  in  depth  since  the  rise  and  fall  rates  are  ~  1  m/sec.  CTD 
data  are  collected  once  per  second,  so  that  they  are  twenty  times  as 
dense  as  the  velocity  data. 

The  horizontal  position  of  the  instrument  is  differenced  and 
provides  a  direct  measurement  of  the  horizontal  velocity  components. 

A  thorough  analysis  of  the  sources  and  magnitudes  of  errors  inherent  in 
or  incidental  to  the  White  Horse  System  (i.e.,  unavoidable  and  random 
errors)  was  performed  by  Luyten,  Needell  and  Thomson  (1982).  A 
hydrodynamical  analysis  of  instrument  response  time  was  dene  by 
Hendricks  and  Rodenbusch  (1981). 

On  occasion,  the  third  bottom  transponder  has  been  recovered  (for 
use  at  another  station)  before  all  instrument  drops  at  a  location  were 
completed.  In  those  cases,  the  triangulation  calculation  proceeds 
using  two  transponders  and  depth,  the  latter  obtained  from  the  White 
Horse  CTD  pressure  sensor.  Absolute  velocity  errors  are  somewhat 
higher,  for  that  case  (cf.  Figure  7  of  Luyten,  Needell  and  Thomson, 
1982).  Table  1 1—2  indicates  which  triangulation  method  was  used  for 
each  profile. 

A  typical  survey  pattern  is  shown  in  Figure  1 1-2.  The  surface 
track  which  the  ship  steamed  while  performing  the  survey  is  indicated 
by  a  dashed  line.  Individual  survey  points  (shown  as  black  dots)  were 
collected  while  the  ship  was  freely  drifting  and  are  used  in  the 
least-squares  fit  to  locate  the  beacons  with  respect  to  each  other. 


1  he  bottom  positions  of  the  transponders,  determined  from  the  survey 
frocessing,  are  shown  as  crossed  circles. 

Figure  1 1—3  presents  profile  6(  7  in  plan  view.  Only  the  horizontal 
location  of  the  White  Horse,  from  release  to  retrieval,  is  shown.  The 
symbol  "+"  marks  the  down  track,  "a"  the  up  track.  The  spacing  of  the 
symbols  indicates  that  the  fall  rate  was  slightly  slower  (about  20  per 
cent)  than  the  rise  rate  so  that  the  instrument  had  less  time  to  move 
horizontally  within  the  flow  field  vhile  rising  vertically.  All  the 
data  used  here  have  been  interpolated  onto  a  uniform  25m  depth  grid,  to 
provide  regular  series  for  comparison. 

Because  the  conductivity  sensor  is  located  at  the  bottom  of  the 
instrument  when  it  rides  in  the  water,  the  only  useful  CTD  data  are 
from  the  down  profiles  when  the  sensor  is  not  affected  by  instrumental 
entrainment.  Vertical  isopycnal  displacements  calculated  from  the  CTD 
data  are  thus  available  as  profiles  corresponding  to  the  down  velocity 
profiles.  Consequently  only  the  down  portions  of  drops  were  used  for 
cross-spectral  estimates  and  for  the  dropped  lagged  coherence 
estimates.  The  exception  is  profile  705,  for  which  only  an  incomplete 
down  trace  was  obtainable.  The  velocity  data  from  up  and  down  profiles 
are  equally  accurate,  however,  and  both  profiles  were  used  for  the 
autospectral  estimates  in  Chapters  III  and  IV.  Comparison  of  up  and 
down  profiles  indicates  that  there  are  no  nontidal  systematic 
differences  (Luyten,  Needell  and  Thomson,  1982). 

Table  1 1-2  summarizes  pertinent  information  about  each  White  Horse 
profile.  In  a  few  indicated  cases,  a  physically  different  but 


POSITION.  METERS 


DROP  «  607  ...LOCATION  OF  WHITEHORSE 


01STRNCE,  IN  METERS.  RELRTIVE  TO  8ERCQN  RT  (0,0) 


PLU5  SYMBOL:-  GOING  OOWN 

SMALL  TMIflNGLE  SYMBOL:-  COMING  UP 


FIGURE  1 1 -3 

White  Horse  profile  track  in  plan  view,  showing  zonal  and 
latitudinal  peregrinations. 
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TABLE  I 1-2 

1976  WHITi:  HORSE  DATA 


Profile 

Date 

and 

Time 

Code 

Beacon 

Notes 

501 

May 

14 

0820 

DSX 

23 

Top  450m  of  up  profile  combined 
with  rest  of  down  profile. 

Signal  lost  below  4765m.  Beacon 
1  unusable  due  to  distance 

502 

May 

28 

1035 

23 

503 

May 

28 

1706 

23 

504 

June 

18 

1010 

BDS 

23 

CTD  flooded  at  2464m  down 

601 

May 

16 

0900 

Z 

23 

Beacon  1  noisy  or  no  signal. 

602 

May 

29 

1238 

* 

1? 

North  comp,  noisy,  upper  500m  up 

603 

June 

4 

1854 

12 

No  data 

604 

June 

5 

0018 

*S 

1? 

605 

June 

6 

0635 

* 

12 

606 

June 

11 

0540 

*D 

12 

Beacon  2  released  from  bottom 

with  dropsonde 

new  net,  due  to  beacon  replacement 


607 

June 

12 

2112 

* 

11! 

608 

June 

15 

1820 

* 

12 

40-second  sampling  interval 
interpolated  to  20  sec. 

609 

June 

15 

1933 

B 

12 

609.. heavy  interpolation, 
less  on  up 

610 

June 

16 

1855 

B 

12 

Heavy  interpolation  on  up,  above 
2490m 

611 

June 

17 

1216 

B 

12 

701 

May 

17 

1321 

A 

Instrument  close  to  baseline 

702 

May 

30 

0815 

* 

23 

Beacon  1  noisy 

703 

June 

4 

0655 

★ 

A 

Instrument  close  to  baseline 

704 

June 

6 

1659 

* 

A 

Instrument  close  to  baseline 

705 

June 

10 

2032 

*4 

A 

Down  starts  at  2775m,  (4) 

Up  noisy,  but  O.K. 

706 

June 

13 

0637 

*D4 

A 

No  down  profile  above  150m 

Up  not  usable  (4) 

Code  Beacon  Configuration 

*  Good  CTD  data  A  All  3  beacons  used  for  position 

8  Black  Horse  instrument  used  calculations. 

D  Only  'down'  profile  available  12  Beacons  1,2  and  pressure.... 

S  Profile  shallower  than  4900  m  23  Beacons  2,3  and  pressure.... 

X  Beacon  1  retrieved  after  drop 
Z  Beacon  3  retrieved  after  drop 

4  4ms  problem,  see  data  report  (Spencer,  O'Neill  and  Luyten,  1980) 


■ 
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TABLE  1 1-2  (Continued) 
1976  WHITE  HORSE  DATA 


Profile 

Date 

and 

1  ime 

Code 

Beacon 

707 

June 

15 

0803 

★ 

A 

708 

June 

17 

0317 

BD 

12 

801 

May 

18 

0254 

A 

802 

May 

30 

2346 

*Z 

A 

803 

June 

3 

1855 

* 

12 

804 

June 

7 

0349 

*D 

12 

805 

June 

10 

0753 

★ 

12 

806 

June 

13 

1517 

★ 

12 

807 

June 

14 

2343 

*D4 

12 

901 

June 

1 

0106 

*Z 

A 

902 

June 

3 

0444 

* 

12 

903 

June 

7 

1818 

*D4 

12 

904 

June 

9 

1821 

* 

12 

905 

June 

14 

0411 

* 

12 

1001 

June 

2 

0643 

*Z 

A 

1002 

June 

8 

1436 

* 

12 

1003 

June 

9 

0125 

8 

12 

1101 

June 

5 

2202 

* 

AC 

1102 

June 

12 

0440 

★ 

A 

1103 

June 

16 

0547 

B 

AC 

Notes 


Heavy  interpolation  above  1600m 
Pressure  case  flooded  at  3200m  up 
No  signal  from  beacon  3 


Up  not  usable,  multiple  problems 
No  down  profile  above  350m 
Noisy  data  upper  200m 
No  down  profile  above  200m 
5  point  smoothing  done  on  noisy 
pressure  data.  Heavy  interpola¬ 
tion  around  750m  down  (4) 

Up  not  usable  (4) 

No  up  profile  above  300m  (rainl) 

No  down  profile  above  100m 
Up  not  usable  (4) 


Heavy  interpolation  300-450m  up 
No  up  profile  above  475m 


Below  3000m,  profiles  formed  from 
2,3  beacon  case 

950-3100m,  profiles  formed  from 
2,3  beacon  ca;e 


Code 

* 

B 

0 

S 

X 

Z 

4 


Beacon  Configuration 

Good  CTD  data  A  All  3  beacons  used  for  position 

Black  Horse  instrument  used  calculations. 

Only  'down'  profile  available  12  Beacons  1,2  and  pressure - 

Profile  shallower  than  4900  m  23  Beacons  2,3  and  pressure.... 

Beacon  1  retrieved  after  drop 
Beacon  3  retrieved  after  drop 

4ms  problem,  see  data  report  (Spencer,  O'Neill  and  Luyten,  1980) 
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electronical 1y  identical  version  of  the  instrument,  n  imed  the  Black 
Horse,  was  used  in  place  of  the  White  Horse.  A  data  -eport  covering 
both  the  White  Horse  profiles  analyzed  here  and  assoc  ated  current 
meter  observations  was  prepared  in  1980,  and  contains  additional 
details  regarding  both  (Spencer,  O'Neill  and  Luyten,  980).  Comparison 
of  White  Horse  and  simultaneous  current  meter  velocity  measurements  has 
teen  used  to  confirm  the  velocity  error  estimates  above  (Luyten, 
taedell  and  Thomson,  1982)  and  also  the  absolute  orientations  of  the 
transponder  nets  (Needell,  personal  communication). 

The  White  Horse  system  has  been  used  in  several  other  locations  and 
c :eanographic  experiments: 

-  Western  Sargasso  Sea  [M0DE-0]  1972  (Gould,  Schmitz  and 

Wunsch,  1974) 

-  Gulf  Stream  1974  (Spencer,  1979,  WHOI  data  report) 

-  Equatorial  western  Pacific  1978  (Eriksen,  1981) 

-  Equatorial  western  Indian  Ocean  1979  [INDEX].  Some  of  the 

same  plus  additional  stations;  much  more  extensive 
coverage  along  the  equator,  from  47°E  to  62°E. 

(Luyten  et  al.,  1980;  Levy,  Spencer  and  Needell,  1981, 
WHOI  data  report,  in  preparation) 

-  Equatorial  central  Pacif  c  1981  [PEQUOD], 

Section  B.  Calculations  from  the  C TP  Data 

B.l.  Buoyancy  Frequency.  Buoyancy  frequency  profiles  were 
calculated  using  a  least  squares  fit  of  the  gradient  of  potential 
specific  volume  anomaly,  for  all  drops  where  CTD  data  were  available 
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and  could  be  edited.  Table  1 1—2  indicates  the  profiles  for  which  CTD 
data  were  available.  [Those  from  Black  Horse  drops  (excluded  because 
its  CTD  malfunctioned),  drops  which  did  not  reach  the  bottom  and  those 
immediately  after  replacement  of  a  bro  en  conductivity  cell  were  not 
used.  There  is  no  CTD  data  from  0°,  5  )° 30 * E ;  all  data  is  from  the  53*E 
meridian.]  An  average  buoyancy  frequency  profile  [=N(x,z)=N(z)]  was 
produced,  in  order  to  minimize  tie  perturbations  on  the  background 
stratification  peculiar  to  a  particular  place  or  time.  This  average 
was  used  to  reduce  each  velocity  profile  to  the  form  which  might  have 
been  observed  in  a  uniformly  stratified  ocean. 

The  temporal  averages  of  N2(x,y,z)  at  each  latitude  were 
ensemble-averaged  to  produce  a  profile  of  x , z) .  The  only  spatial 
smoothing  applied  is  thus  latitudinal.  No  smoothing  was  done  over 
depth.  Such  smoothing  was  examined  and  had  no  visual  or  spectral 
impact.  The  buoyancy  frequency  prof i If  used  for  WKB  stretching  and 
normalizing  as  well  as  for  the  displacunent  calculations  (below)  is 
shown  in  Figure  I 1—4. 

B.2.  Vertical  Displacement.  Vertical  displacements  of  isopycnals 
from  their  equilibrium  positions  have  teen  calculated  by  a  comparison 
of  individual  CTD  profiles  to  a  space-time  average  of  sigma  theta  which 
defines  the  background  stratification.  The  linearized  conservation  of 
mass  equation  for  an  incompressible  Boussinesq  fluid  can  be  written 


3p 


0 


3t 


3P, 


az 


(II-l) 


where  p  is  the  deviation  of  potential  density  from  the  mean  (p  -  p  -p  ; 

v  V  v  O 


FIGURE  I 1—4 


Plot  of  N(z)  against  depth:  illustraticn  of  the  average 
buoyancy  frequency  profile  used  for  WKB  stretching  and 
normalizing,  and  in  the  calculation  of  vertical 
displacement. 
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p0  is  the  observed  potential  density  [  =  1  +  10"^oe]).  Since  ]  =  J wdt, 
where  ^  is  verticil  displacement  measured  upward,  equation  Il-l  can  be 
integrated  in  time  to  yield 


]  (x,y,z,t) 


Pe(x,y,z,t)  -  p9(x,z) 
- - 


__a_ 

N2(z) 


10 


( 1 1-2) 


~7  q  pe 

where  N  (z)  =  — *  — —  .  Equation  1 1—2  was  used  to  calculate  the 

o  Z 

pe 

/v 

displacement  profiles  presented  here.  The  computation  of  oa  included 
a  least  squares  fit  of  the  in  situ  potential  specific  volume  anomaly 
over  a  variable  depth  interval,  referred  to  the  center  of  the 
interval:  lOOdb  above  500m,  calculated  every  25db;  2Q0db  between  500 
and  4000m,  and  400db  below  4000m,  calculated  every  50db  for  the  latter 
intervals.  Values  above  50m  and  below  -4800m  are  extrapolated,  and 
the  profile  is  interpolated  to  the  same  25m  depth  grid  used  for  the 
velocity  data.  Adjacent  values  are  therefore  not  independent  anywhere 
in  the  water  column,  which  accounts  for  the  relative  smoothness  with 
depth  of  the  variability  in  the  displacement  profiles.  (See  Appendix 
A.)  Average  was  obtained  by  ensemble-averaging  the  temporal 
averages  at  latitudes  3/4°S  to  3*N.  The  buoyancy  frequency  profile 
used  to  compute  displacement  is  the  space-time  average  discussed  above 
and  used  below  for  WKB  stretching  and  normalizing, 
he  c<  lculation  can  be  visualized  graphically: 


■<*  *’ 
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S  .etch  illustrating  calculation  of  vertical  isopycnal 
d  splacement  using  deviations  of  an  individual  sigma 
tueta  profile  from  the  space-time  average  sigma  theta 
profile. 


.  Vzl)  ' 

i.  ](.-,)  =  - = - —  .  Density  fluctuations  due  to  horizontal 

<V21} 

~ Tz - 

di ‘.pi aceme.its  have  been  ignored,  and  the  underlying  assumption  is  that 
spatial  av  raging  over  four  degrees  of  latitude  and  temporal  averaging 
over  three  weeks  provides  a  sufficient  approximation  to  the  equilibrium 
density  fit  Id.  For  sample  profiles,  see  Figure  1 1—6.  The  full  set  of 
displacemeit  profiles  are  shown  as  latitudinal  sections  in  Appendix  A 
an  I  as  temporal  sections  at  each  station  in  Appendix  B. 

Section  C.  WKB  Stretching  and  Normalizing  Procedures 

To  permit  simultaneous  treatment  of  upper  and  lower  portions  of  the 
water  column,  the  apparent  dependence  of  velocity  amplitude  and 
vertical  length  scale  upon  the  buoyancy  frequency  has  been  removed, 
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FIGURE  II -6 

Sample  of  displacement  profiles.  These  are  the  five 
profiles  available  at  the  equator,  53*E.  See  Appendices 
A  and  B  for  full  set  of  profiles. 
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using  a  WKB  stretching  and  normalizing  procedure  first  discussed  by 
Bell  (1974).  Observational  evidence  justifying  this  procedure,  based 
on  moored  current  meter  measurements,  can  be  found  in  Luyten  (1982). 
Theoretical  justification  is  based  upon  the  solutions  to  the  vertical 
structure  equation  discussed  in  Chapter  I-B. 

Horizontal  velocities  are  scaled  by  dividing  them  by  [N(z)/N0]l/2, 
where  N(z)  is  the  average  buoyancy  frequency  profile  described  in 
Section  B.l  and  NQ  is  the  average  of  that  profile.  The  NQ  used  here 
is  1.0915  cph,  for  reasons  which  are  mentioned  below. 


The  vertical  coordinate  is  correspondingly  stretched,  by  multiplying 
each  depth  increment  by  the  buoyancy  frequency  at  that  level,  normalized 

hy  n0. 


When  the  NQ  used  is  the  average  of  the  Brunt-Vaisala  profile  used,  as 
here,  the  stretched  vertical  profile  (u*,  v*  or  J*  as  functions  of  z*) 
maps  back  onto  the  original  total  water  depth,  5000m  in  this  case 
becoming  5000  stretched  meters  (sm).  Vertical  wavelengths  and  scales 
can  then  be  readily  interpreted  as  fractions  of  the  total  water  depth, 
i.e.,  xz«1200sm  implies  approximately  four  wavelengths  between  the 
surface  and  the  bottom. 
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The  maximum  stretched  depth  interval  between  25m-data  points  is 
263sm  at  100m,  where  N  is  13.4  cph.  The  maximum  interval  in  the  depth 
series  used  for  spectral  analysis  (see  Section  D)  is  63sm,  corres¬ 
ponding  to  N(262.5m)=2. 7  cph.  The  minimum  stretched  depth  interval  is 
4sm,  where  N(4987.5m)  is  0.18  cph.  The  unequally-spaced  stretched  data 
have  been  interpolated  onto  a  uniform  25sm  depth  grid. 

An  adaptive  filtering  technique  was  used  to  smooth  the  stretched 
data  in  the  lower  part  of  the  water  column,  where  a  100m  interval  in 
the  vertical  coordinate  is  condensed  into  as  little  as  16sm.  A  boxcar 
filter  was  applied  to  the  scaled  data  before  stretching.  The  filter 
length  is  determined  by  the  relationship  between  the  local  and  the 
depth-averaged  buoyancy  frequencies,  according  to  the  following  formula 
(Dunlap,  personal  communication): 

z  *  naz 


x*(z)  = 


2n 


~h  £ 


U)  , 


(1 1-3) 


z  -  n&z 

where  n  =  1/2  NQ/N(z).  The  filter  has  no  impact  upon  the  upper  part 
of  the  stretched  water  column,  where  the  local  buoyancy  frequency  is 
greater  than  NQ,  but  smooths  data  in  the  lower  portion  over  as  much 
as  150sm. 

The  resultant  velocity  vectors  have  been  used  for  the  analysis  to 
be  described.  The  upper  300m  of  the  water  column  is  expanded  to  the 
upper  1500sm  of  the  stretched  profile;  the  bottom  kilometer  of  the 
water  column  is  condensed  to  the  lower  300sm  of  the  stretched  profile. 
Upper  level  velocities  have  been  relatively  reduced  and  lower  level 
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velocities  enhanced.  One  can  now  discuss  vertical  wavenumbers 
throughout  the  temporal  and  spatial  array  rather  than  merely  local 
vertical  wavenumbers.  Figure  I I— 7  shows  the  stretched  vertical 
coordinate  (z*)  plotted  against  depth  (z)  and  illustrates  the  one-to- 
one  mapping  of  5000m  onto  5000sm.  The  upper  part  of  the  figure, 
corresponding  to  the  top  275m  of  the  total  ocean  depth,  has  been 
crosshatched,  to  indicate  the  equivalent  top  I400sm  of  the  stretched 
coordinate.  The  impact  on  the  original  data  of  the  WKB  stretching  and 
normalizing  is  evident  in  Figure  1 1—8. 

Section  D.  Spectral  Analysis  Methods 

Time  series  analysis  has  been  used  in  this  study  in  two  ways: 
conventionally,  as  a  means  of  determining  the  distribution  of  energy 
among  vertical  wavenumbers  (autospectra)  and  phase  relationships 
between  variables  (cross-spectra),  and  less  conventionally,  looking  at 
dropped  lagged  coherence  and  phase  (a  particular  averaging  of 
cross-spectra)  to  detect  vertically  standing  or  propagating  waves.  In 
both  cases,  the  upper  275m  of  the  water  column  has  been  ignored.  The 
density  gradient  is  changing  rapidly  there  and  the  first-order  WKB 
ipproximation  implicit  in  the  scaling  is  not  applicable.  In  addition 
(and  perhaps  more  importantly),  linear  wave  theory  is  unlikely  to  be  an 
appropriate  model  in  that  region,  so  that  spectral  decomposition  of 
this  portion  of  the  data  would  be  less  informative.  The  total  depth 
interval  analyzed  is  3600sm,  so  that  the  maximum  spectral  resolution 
with  standard  techniques  is  wavelengths  of  3600,  1800,  1200,  900, 

...sm.  The  top  of  the  series  used  is  1412. 5sm,  correspond’’ ng  to  the 


r  i  -vi  r--  - 


FIGURE  I I -7 


Plot  of  2*  against  z:  stretched  vertical  coordinate 
versus  depth.  Illustrating  the  smoothness  of  the  WKB 
one-to-one  mapping. 
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original  data  point  262. 5m( 1399sm).  The  deepest  point  used  is 
4987. 5sm,  corresponding  to  4887.5m(4990sm). 

A  10%  1/2-cosine  bell  is  used  to  taper  the  data  prior  to  spectral 
decomposition  (Koopmans,  1974,  p.302).  This  taper,  a  smooth  modulating 
function,  removes  end  effects,  but  also  removes  half  the  energy  from 
the  first  and  last  ten  percent  of  the  record.  The  final  spectral 
values  should  thus  be  adjusted  by  a  factor  of  1.11  (Munk,  Snodgrass  and 
Tucker,  1959),  before  comparison  with  untapered  spectra. 

D.l.  Auto-  and  cross-spectra.  A  question  which  can  be  addressed 
via  autospectral  analysis  is  that  of  vertical  homogeneity.  Decom¬ 
position  of  the  3600sm  profiles  at  net  6  into  upper  and  lower  1800sm 
pieces  reveals  no  significant  differences  in  their  autospectra.  The 
two  halves  were  not  tapered  prior  to  the  spectral  computation. 

Statistical  reliability  of  spectral  estimates  with  no  wavenumber 
averaging  is  dependent  upon  the  number  of  pieces  (in  this  case, 
profiles  available  at  each  latitude)  included  in  each  estimate.  CTD 
data  are  available  for  a  subset  of  profiles  at  each  net,  and  have  been 
used  to  compute  vertical  displacement  profiles.  At  the  equator,  CTD 
data  and  thus  vertical  displacement  are  available  for  only  half  of  the 
velocity  profiles,  spanning  a  total  of  two  and  a  half  weeks.  At  3/4*N 
CTD  data  are  available  for  six  of  the  eight  drois,  spanning  the  same 
time  period  as  at  the  equator.  The  maximum  numuer  of  vertical 
displacement  profiles  at  a  single  latitude  is  six  at  3/4 "N  and  1.5*N, 
and  falls  to  a  minimum  of  two,  at  3/4*S  and  5*N.  There  are  five 
profiles  at  two  stations,  the  equator  and  3*N.  All  of  the  above  are 
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along  the  53*E  meridian.  There  is  none  available  at  the  other 
equatorial  station  at  50*30’E. 

A  relevant  question  regarding  the  cross-spectral  phase  estimates 
involves  the  accuracy  of  the  determination  of  the  absolute  net 
orientations.  It  was  asserted  in  Section  A.l  that  the  orientations  had 
been  resolved  to  within  ±  5*.  General  conclusions  from  a  test  of  the 
impact  of  orientation  errors,  by  rotation  of  the  velocity  vectors,  are 
that  the  combination  of  high  coherence  and  a  90*  phase  estimate  is 
invariant  with  respect  to  rotation,  since  that  result  is  dominated  by 
the  quadrature  spectrum,  which  is  invariant  (Fofonoff,  1969).  The 
combination  of  high  coherence  and  a  ~180*  phase  estimate  js^  affected  by 
rotation,  however.  Rotation  in  either  direction  produces  at  first  a 
gradual  decrease  in  both  coherence  and  phase.  A  rotation  of  30°,  for 
instance,  changes  the  phase  estimate  for  net  6  at  a  wavelength  of 
1200sm  by  only  16*,  while  the  coherence  remains  significantly  non¬ 
zero.  However,  coherence  amplitude  drops  to  a  minimum  at  45*,  when  the 
components  are  in  quadrature,  and  then  increases  to  a  maximum  at  90*, 
when  the  components  are  in  phase.  It  seems  unlikely,  therefore,  that 
any  of  the  cross-spectral  results  presented  in  Chapters  III  and  IV  are 
affected  by  net  orientation  errors. 

D.2.  Dropped  Lagged  Coherence  and  its  interpretation.  Since  the 
1976  White  Horse  velocity  profiles  constitute  a  time  series  of 
snapshots  of  the  water  column,  it  is  possible  to  use  them  to  examine 
temporal  behavior  as  well  as  vertical  scales.  Specifically,  it  is 
possible  to  use  them  to  examine  the  question  of  vertical  phase 
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propagation  of  the  observed  oscillations.  The  method  applied  is  known 
as  dropped  lagged  coherence  (OLC)  (cf.  Hayes,  1975).  Pairs  of  profiles 
separated  by  the  same  temporal  lag  are  grouped  and  their  ensemble- 
averaged  cross-spectra  computed.  For  each  variable,  changes  in 
coherence  and  phase  estimates  with  temporal  lag  can  then  be  analyzed  in 
a  particular  vertical  wavenumber  band,  and,  in  special  cases,  lend 
themselves  to  straightforward  interpretation. 

An  assumption  underlying  the  approach  of  dropped  lagged  coherence 
is  that  vertical  profiles  at  the  same  geographical  location  represent 
stationary  processes,  so  that  the  cross-spectra  for  a  particular 
wavenumber  band  are  a  function  primarily  of  the  time  difference  between 
the  profiles.  The  cross-spectral  estimates  for  each  vertical  wave- 
number  and  for  each  pair  of  profiles  are  assumed  to  be  independent.  A 
certain  degree  of  homogeneity  within  a  profile  is  required  in  order  to 
produce  a  detectable  signal,  i.e.,  a  single  process  mist  be  responsible 
for  variability  throughout  a  substantial  portion  of  the  water  column  at 
a  particular  vertical  scale. 

Two  simple  models  of  wavelike  motion  in  vertical  profiles  are 
Illustrated  in  Figure  1 1—9.  For  particular  frequency  and  vertical 
wavenumber,  the  models  are  of  a  standing  wave  (vertical  mode)  and  of  a 
propagating  wave.  For  a  standing  wave,  dropped  lagged  coherence  would 
drop  to  zero  and  phase  would  change  by  ir  at  temporal  lags  equal  to 
1/4  and  3/4  of  the  wave  period.  Coherence  would  be  high  at  temporal 
lags  of  some  integral  multiple  of  half  the  wave  period.  For  a 
vertically  propagating  wave,  coherence  would  be  expected  to  maintain  a 
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constant  value  (probably  less  than  unity  because  of  noise,  and  possibly 
diminishing  with  temporal  lag  due  to  processes  such  as  dissipation,  for 
example),  and  phase  would  be  proportional  to  the  time  lag,  progressing 
linearly  through  a  full  360*  as  the  temporal  lag  approached  the  period 
length.  (4  *  «t,  where  4  is  the  phase  estimate  and  «  is  the  wave 
frequency.)  The  sign  of  the  phase  relationship  would  then  indicate  the 
direction  of  vertical  phase  propagation. 

The  absence  of  coherence  or  the  presence  of  random  phase  relations, 
on  the  other  hand,  would  indicate  that  neither  of  these  simple  models 
was  appropriate  or  demonstrable. 

Only  a  restricted  range  of  periods  can  be  inferred  from  DLC  of  the 
1976  White  Horse  profiles.  There  are  eight  resolvable  lags  at  the 
equator,  to  a  maximum  lag  of  19  days,  and  only  six  at  3/4°N,  out  to  17 
days.  The  gaps  between  the  first  and  second  profiles  at  each  latitude 
are  too  large  to  allow  an  examination  of  longer  lags.  Dropped  lagged 
coherence  was  computable  at  1.5*N  as  well,  although  only  iour  temporal 
lags  are  resolvable.  The  total  time  span  of  the  series,  one  month,  is 
too  short  to  resolve  reliably  periods  longer  than  several  months.  The 
total  number  of  profiles  at  one  location  and  their  spacing  in  time 
constrains  both  the  number  of  lags  resolvable  and  the  gaps  between  them 
and  affects  the  range  of  elapsed  time  collected  at  a  single  lag, 
putting  a  lower  limit  of  several  weeks  on  the  periods  that  might  be 
observed. 

In  the  context  of  linear  equatorial  wave  theory,  the  period  range 
of  several  weeks  to  several  months  resolvable  using  the  dropped  lagged 
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coherence  estimates  means  that  only  Kelvin  and  mixed  Rossby-gravity 
waves  are  likely  to  be  detectable  using  the  method  on  this  depth 
series.  Table  1-2  showed  the  maximum  period  at  which  inertia-gravity 
waves  can  exist  for  each  of  the  wavelengths  resolvable,  and  the  minimum 
period  at  which  Rossby  modes  can  exist.  Only  Kelvin  and  mixed 
Rossby-gravity  waves  satisfy  the  linear  dispersion  relation  in  the 
period  range  between  them:  between  15  and  86  days,  for  example,  for  a 
vertical  wavelength  of  1200sm. 
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CHAPTER  III  -  Zonal  Velocity 


Introduction 

The  dominant  characteristics  of  the  zonal  velocity  field  from  these 
profiles  are  the  persistence  with  time  and  latitude  of  regular  features 
in  the  vertical  structure,  and  the  phase  reversal  observed  between 
profiles  at  1.5°N  and  3*N.  The  approach  in  this  chapter  will  be  to 
describe  and  quantify  the  distribution  of  zonal  kinetic  energy  with 
latitude  (Section  A),  to  examine  relationships  between  zonal  velocity 
and  other  variables  (Section  B),  to  examine  evidence  of  periodicity  in 
Section  C,  and  to  evaluate  possible  interpretations  in  terms  of  linear 
iquatorial  wave  theory  in  Section  D.  Comparison  with  other  results 
will  be  made  in  Chapter  V. 

Section  A.  Latitudinal  Distribution  of  Zonal  Velocity 

A. 1.  Latitudinal  Energy  Density,  The  latitudinal  distribution  of 
zonal  kinetic  energy  as  a  function  of  vertical  wavenumber  is  shown  in 
Figure  1 1 1— la.  four  of  the  curves  are  labelled  by  the  center 
wavelength  (in  stretched  meters)  of  the  averaging  interval:  each  point 
plotted  represents  an  average  over  five  wavenumber  bands.  The 
unlabelled  curve  at  the  top  of  the  figure  is  the  low-order  estimate  and 
represents  an  average  of  all  wavelengths  greater  than  1800sm,  i.e.,  the 
profile  mean,  all  trends,  and  wavelengths  as  short  as  half  the  examined 
water  column.  The  error  bars  on  the  estimates  comprising  the  labelled 
curves  represent  the  95%  confidence  interval  limits. 

A  concentration  of  energy  near  the  equator  is  evident  in  the 
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shortest  wavelength  bands  and  is  particularly  strong  for  2 40sm.  The 
energy  density  estimates  at  the  equator  and  1.5°N  are  significantly 
different  from  each  other  for  all  four  wavelength  bands  presented  here, 
from  720sm  through  180sm.  Also  notice  that  the  estimates  for  the  720sm 
and  the  360sm  band  at  each  latitude  are  separated  by  an  order  of 
magnitude  in  energy  density.  The  estimates  for  360,  240  and  180sm 
cluster  one  and  a  half  orders  of  magnitude  below  that  for  720sm. 
Comparison  of  Figure  1 1 1— la  with  the  corresponding  figure  for  the 
meridional  velocity  component  (Figure  I V— 1 )  reveals  that  both  the 
estimated  energy  levels  and  the  concentration  of  energy  density  at  the 
equator  for  the  horizontal  velocity  components  are  not  significantly 
different  from  each  other  for  the  larger  wavenumbers  examined  here: 

360,  240  and  180sm.  However,  the  shapes  of  the  curves  with  latitude 
and  their  levels  are  significantly  different  (at  the  95^  confidence 
level)  for  the  two  shortest  wavenumber  bands,  those  for  a  centered 
wavelength  of  720sm  and  for  the  low-order  estimate  (uppermost  curve). 

Figure  1 1 1— lb  shows  the  similarly  averaged  autospectral  estimates 
for  the  vertical  displacement  field,  by  latitude.  The  pressure  wave 
field  for  equatorial  waves,  and  thus  the  displacement  field,  has  no 
simple  structure  about  the  equator,  such  as  the  Gaussian  velocity 
structures  associated  with  the  Kelvin  and  mixed  Rossby-gravity  waves. 
The  patterns  here  certainly  do  not  indicate  equatorial  trapping.  The 
estimates  for  a  particular  wavenumber  band  are  indistinguishable  at  the 

a 

95  confidence  level. 


Autospectral  estimates, 


FIGURE  Ill-lb 

t>y  latitude,  for  vertical  displacement. 
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A. 2.  Latitudinal  Profile  Sections.  Figure  1 1 1—2  illustrates  both 
the  persistence  of  vertical  structure  features  across  the  equator  and 
as  far  off  the  equator  as  1.5‘N,  and  also  the  striking  antisymmetry 
between  the  vertical  structure  of  zonal  velocity  at  1.5°  and  3*N.  The 
profile  data  (u,  u*,  ]  ,  and  1*)  are  taken  from  latitudinal 
section  =*3,  spanning  June  2  to  June  5,  1976.  The  complete  set  of 
latitudinal  sections  is  presented  in  Appendix  A.  The  equatorial 
vertical  displacement  profile  is  not  part  of  the  displacement  subset 
jsed  here,  since  for  this  meridional  section  the  CTD  data  at  the 
equator  was  incomplete. 

A  simple  measure  of  the  degree  of  symmetry /anti symmetry  between 
.hese  profiles  is 


where  Uj  and  U2  are  the  zonal  velocities  of  two  profiles  being 

compared.  For  a  case  of  antisymmetry.  Equation  III-l  will  yield  small 

values  of  P  as  one  profile  reflects  into  the  other:  the  residual 

variance  is  then  the  variance  which  cannot  be  explained  by  pure 

reflection.  In  the  opposite  case  of  strong  symmetry,  P  will  be  close 

to  2.0.  A  value  of  one  for  P  is  thus  analogous  to  a  cross-correlation 

coefficien.  of  zero,  indicating  that  neither  symmetry  nor  antisymmetry 

cominates,  and  correspondingly,  P  =  2  is  similar  to  R..  =  1; 

ulu2 

P  =  0  R  =  -1. 

U1UL 

Computing  P  with  the  profiles  from  3#N  and  1.5#N  in  this  manner 
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(summing  over  the  bottom  3600sm  of  the  stretched  water  column),  gives 
the  result  that  almost  60%  of  the  variance  can  be  removed  by 
reflection.  Neglect  of  the  top  1400sm  was  discussed  in  Chapter  1 1— D  in 
relation  to  the  applicability  of  the  WKB  approximation  implicit  in  the 
stretching  and  normalizing  procedure  applied  to  the  data.  [When  the 
entire  stretched  water  column  is  employed  in  the  computation  of  P  , 
the  result  is  that  only  24%  of  the  variance  is  removed.  However,  as 
discussed  in  Chapter  II,  the  upper  ocean  currents  are  not  well-resolved 
in  the  stretching  and  normalizing  procedure:  the  upper  300m  of  the 
water  column  has  been  expanded  to  the  upper  I500sm  of  the  stretched 
water  column.  In  addition,  this  portion  of  the  water  column  is 
expected  to  be  nonlinear.] 

Table  III-l  present:  the  results  of  the  computation  of  Equation 
III-l  for  various  profile  pairs  from  meridional  section  #  3.  The  upper 
portion  of  the  table  compares  profiles  whose  antisymmetry  can  be  seen 
visually  in  Figure  1 1 1—2.  In  general,  38-59%  of  the  variance  in  pairs 
of  profiles  spanning  2*N  can  be  explained  in  terms  of  antisymmetry. 

The  lower  portion  of  the  table  indicates  the  degree  of  persistence  with 
latitude  of  structural  features  in  the  zonal  velocity  profiles.  From 
40-63%  of  the  total  variance  is  removed  via  subtraction  of  one  profile 
from  another.  Note  that  the  least  variance  is  removed  when  the 
equatorial  profile  is  compared  to  a  profile  one  and  a  half  degrees 
away,  a  fact  related  to  the  shorter  vertical  scales  present  at  and  near 
the  equator  (viz.  Figure  III-l  and  the  upper  part  of  Figure  1 1 1—2) . 

Also  note  that  the  last  line  of  the  table  examines  the  degree  of 
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TABLE  1 1 1 — 1 :  Quantif icatioi  of  symmetry/antisymmetry  of  zonal 
velocity  profiles  with  latitude  (using  bottom 
3600sm  of  each  profile  unless  otherwise  noted). 
All  profiles  are  from  latitudinal  transect  3. 
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.29 
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.59 
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symmetry  across  the  equator,  by  comparing  the  profile  at  3/4°S  with 
that  at  3/4*N.  Symmetry  accounts  for  58%  of  the  variance,  in  this 
case.  The  important  point  here  is  that  neighboring  profiles  that  span 
2*N  are  as  negatively  correlated  as  the  profiles  equatorward  of  2*  are 
positively  correlated. 

A. 3.  Variance -preserving  net  autospectra.  The  autospectral 
estimates  for  zonal  velocity  of  Figure  1 1 1— la  are  presented  in  greater 
vertical  wavenumber  detail  in  Figure  1 1 1—3.  Each  curve  is  for  a 
particular  latitude;  all  profiles  at  each  latitude  along  53°E  were 
averaged  together,  but  no  averaging  was  done  over  wavenumber  bands. 

The  estimates  at  the  equator  are  averages  for  ten  profiles,  and  at 
3/4 *S  and  5*N  for  only  three.  The  format  is  that  known  as 
variance-preserving,  since  the  area  under  the  curve  is  directly  related 
to  the  energy  in  the  wavenumber  band.  The  ordinate  is  energy  density 
times  vertical  wavenumber,  plotted  against  log  wavenumber  along  the 
abscissa.  This  type  of  spectral  representation  compensates  for  the 
increasing  density  of  spectral  estimates  per  unit  log  wavenumber  and  is 
more  accessible  to  visual  inspection.  Note  that  an  m-*  spectrum 
would  appear  as  a  straight  horizontal  line,  so  that  the  structure  in 
these  plots  is  a  measure  of  deviation  from  a  (-l)-slope  on  conventional 
spectral  plots. 

The  strongest  single  signal  is  the  peak,  significant  at  the  95 % 
confidence  level,  for  a  wavelength  of  720sm  at  the  equator  and  3/4#N. 
The  secondary  peak  which  appears  for  360sm  at  those  two  profile 
stations,  but  not  at  1.5*N,  is  one  explanation  for  the  lower 
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FIGURE  I I 1-3 


Autospectral  estimates  of  zonal  velocity  -  no  averaging. 
Each  curve  represents  a  different  latitude. 
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correlation,  discussed  in  the  previous  section,  when  the  profile  at  the 
equator  is  compared  to  profiles  more  than  one  degree  from  the  equator. 
These  peaks  for  360sm  are  significant,  as  is  the  peak  at  the  equator 
for  a  wavelength  of  277sm. 

Since  the  area  under  these  curves  represents  the  energy  in  each 
wavenumber  band,  the  "redness"  of  the  zonal  velocity  speitral  estimates 
is  apparent  in  this  figure,  and  is  prevalent  at  all  latitudes 
examined.  For  contrast,  refer  to  Figure  I V— 3,  where  meridional 
velocity  autospectral  estimates  are  displayed.  The  dominant  feature  in 
that  case  is  the  repetition  of  a  peak  at  1200sm,  a  low  but  not  the 
lowest  wavenumber. 

When  the  3600sm  depth  series  is  decremented  at  the  top  and  bottom, 
to  a  3200sm  series,  the  effect  of  variations  in  the  discretization  of 
vertical  wavelengths  can  be  observed.  The  major  difference  in  the 
autospectra  is  a  drop  in  zonal  kinetic  energy  at  a  vertical  wavelength 
of  1067sm,  relative  to  1200sm.  At  the  equator,  energy  was  evenly 
dstributed  between  1200  and  the  adjacent  bands,  while  energy  at  1067  is 
only  20%  that  in  adjacent  bands.  At  3/4*N,  energy  at  1200  was  30-70% 
that  in  adjacent  bands,  while  energy  at  1067  is  5-20%  that  in 
neighboring  bands.  At  both  latitudes,  the  trough  in  zonal  kinetic 
energy  at  a  vertical  wavelength  of  1067sm  is  statistically  significant 
at  the  95  %  confidence  level. 

Section  B.  Relationships  Between  Zonal  Velocity  and  Other  Variables. 

8.1.  General  cross-spectral  results.  Cross-spectral  estimates  were 
computed  at  each  latitude  for  all  vertical  wavelengths  between  the 


three  variables,  zonal  and  meridional  velocity  and  vertical 
displacement.  All  significant  coherences  are  listed  in  Table  1 1 1-2, 
which  also  appears  as  Table  I V— 1.  Some  of  the  results  are  discussed  in 
Chapter  IV  since  their  implications  are  primarily  regarding  the 
structure  of  the  meridional  velocity  component 

"Significant"  is  used  here  and  in  the  following  to  mean 
significantly  different  from  zero  true  coherence  at  the  95^  confidence 
level.  Error  bars  on  the  phase  estimates  indicate  the  95  %  confidence 
limits.  "Leads"  is  used  in  the  table  and  the  text  following  to  mean 
leads  by  u/2  with  depth  (within  the  error  bars  at  the  95  %  confidence 
level).  "Out  of  phase"  means  a  phase  estimate  of  “180°. 

Cross-spectral  phase  estimates  of  0,  rr/2  or  it  lend  themselves  to 
simple  interpretation  in  terms  of  single  equatorial  meridional  wave 
modes  with  particular  vertical  and  meridional  structure.  Left  open  in 
these  cases  are  the  questions  of  zonal  and  temporal  structure. 

However,  phase  estimates  significantly  different  from  these  key  values 
may  indicate  that  either  a  mix  of  frequencies  is  present,  obscuring 
detection  of  individual  wave  components,  or  that  the  approach  is 
misguided  -  that  the  mix  in  all  dimensions  is  too  complex  for  this 
attempt  at  decomposition. 

The  time  series  of  profiles  at  the  equator  and  the  two  northern 
near-equatorial  stations  are  longer  than  at  the  other  three  stations 
along  53*E,  so  that  more  reliable  statistical  estimates  are  possible 
there.  The  energy  trapping  scales  anticipated  from  linear  equatorial 
wave  theory  lead  one  to  expect  more  significant  coherence  estimates  at 
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the  equator  and  at  3/4°N,  for  the  shorter  vertical  wavelengths. 

Symmetry  relationships  between  variables  mitigate  against  obtaining 
elucidating  cross-spectral  estimates  from  the  equatorial  data  for 
linear  wave  modes,  however  (see  Chapters  I— B. 2.d  and  II-D.l).  Recall 
from  Table  1-1  that  the  equatorial  Rossby  radius  of  deformation  is 
equal  to  3/4°  for  a  vertical  wavelength  of  ~5Q0sm.  Consequently,  the 
latitudinal  structure  of  energy  at  wavelengths  much  shorter  than  that 
is  unlikely  to  be  resolved  using  this  data. 

B. 2.  400sm  cross-spectra.  A  consistent  pattern  emerges  from  the 
cross-spectral  estimates  for  a  vertical  wavelength  of  400sm,  however, 
particularly  at  3/4°N.  At  that  latitude,  u  leads  v  and  J  leads  u, 
both  by  ir/2,  so  that  v  and  5  are  180°  out  of  phase.  These  phase 
relationships  are  consistent  with  the  presence  of  a  mixed 
Rossby-gravity  wave  at  that  vertical  wavelength,  but  the  zonal  kinetic 
energy  in  the  400sm  wavelength  band  is  not  above  the  background  level 

viz.  the  autospectral  estimates  in  Figure  1 1 1—3  and  the  related 
c iscussion) . 

Section  C.  Dropped  Lagged  Coherence  Results. 

C. l.  General.  Figure  III -4  is  a  superposition,  for  the  zonal 
velocity  component,  of  the  dropped  lagged  coherences  and  phases  for 

< ach  temporal  lag  resolvable  at  the  equatorial  station  at  53°E.  The 
cverall  pattern  is  one  of  high  cohe 'ence  and  zero  phase  estimates  at 
the  longer  wavelengths.  For  the  shortest  vertical  wavenumber  bands 
resolvable,  the  dropped  lagged  coherence  calculations  yield  significant 
coherence:  for  wavelengths  of  1800,  900,  and  720sm,  at  least  seven  of 
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FIGURE  II 1-4 

Superposition  of  OLC  for  the  zonal  velocity  component  at  the  equator, 
for  all  temporal  lags  resolvable. 

Key:  dark  to  light  corresponds  to  shorter  to  longer  lags. 
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the  eight  resolvable  temporal  lags  at  the  equator  yield  coherence 
estimates  significant  at  the  95%  confidence  level.  At  3/4*N, 
coherence  between  the  lagged  zonal  velocity  components  is  also  high: 
significant  for  2/3  of  the  six  lags  for  1800sm,  5/6  for  120Gsm,  and  at 
all  lags  for  720sm.  For  none  of  these,  however,  are  the  phase 
estimates  signif icantly  different  from  zero,  which,  as  explained  in 
Chapter  II-D.2,  might  be  due  to  either  of  two  factors:  the  energy 
observed  at  those  vertical  wavelengths  is  related  to  the  presence  of 
vertically  standing  rather  than  propagating  waves,  or  the  energy  is  due 
to  coherent  motion  at  longer  periocs  than  can  be  resolved  with  this 
month-long  set  of  observations. 

C. 2.  400sm  wavelength  band.  The  only  vertical  wavelength  for  which 
:he  DLC  calculations  at  the  equator  yield  both  sufficiently  high 
coherence  to  be  significant  and  phase  estimates  which  are  significantly 
non-zero  is  400sm  (Figure  1 1 1—5) .  The  range  of  periods  which  are 
consistent  with  the  phase  change  with  temporal  lag  is  105-216  days, 
with  phase  propagation  upward.  At  3/4’N,  the  400sm  band  is  coherent 
for  half  the  temporal  lags  examined,  and  two  of  the  phase  estimates  are 
nonzero.  The  phases  are  uniformly  negative,  indicating  phase 
propagation  in  the  opposite  direction  from  that  at  the  equator,  but 
encompassing  the  same  period  range  (47-295  days  in  this  case). 

Section  D.  Interpretation  and  Examination  in  Terms  of  Linear  Theory. 

D. l.  Equatorial  trapping.  The  equatorial  concentration  of  zonal 
kinetic  energy  in  Figure  III-l  can  be  examined  in  terms  of  the  Rossby 
radius  of  deformation  (L^).  If  Kelvin  waves  dominate  the  zonal 
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velocity  field  near  the  equator,  the  latitudinal  distribution  of  energy 
will  scale  with  L^.  Table  I I 1—3  contains  a  comparison  of  the 
energy  ratios  between  3/4°N  and  the  equator  to  the  ratios  consistent 
with  a  deformation  radius  scaling,  for  the  averaged  vertical  wavelength 
bands  centered  at  720,  360,  and  240sm,  where  equatorial  intensification 
appeared  strongest.  The  indication  from  the  table  is  that  Kelvin  waves 
are  not  dominant  at  these  vertical  scales. 


Table  III-3: 

Equatorial  trapping  of  zonal 

velocity 

-R 

(km) 

z 

(sm) 

Autospectral 

O' 

estimates 

3/4‘N 

Energy  ratios 
est.3/4*N  „-,3/4°N 

est.O  lR 

98 

720 

238(180-340) 

236(179-339} 

.99 

.48 

69 

360 

36(28-50) 

33(25-46) 

.92 

.24 

56 

240 

19(14.5-26) 

15.5(12-22) 

.82 

.11 

(  )  : 

indicates  95  confidence  limits. 

D.2.  Evaluation  of  latitudinal  symmetry/ anti  symmetry.  The 
strongest  signals  from  visual  inspection  of  the  1976  White  Horse 
profiles  are  1)  the  temporal  persistence  of  vertical  structures  in  the 
zonal  velocity  component,  at  the  equator  in  particular  (see  Appendix 
B),  and  2)  the  combination  of  symmetry  across  the  equator  and 
antisymmetry  between  the  profiles  at  3*N  and  the  near-equatorial 
profiles,  including  the  profiles  at  1.5*N.  There  were  six  latitudinal 
transects  made  along  53*E  in  1976  which  include  the  station  at  3*N, 
shown  as  Sections  2  through  7  in  Appendix  A.  Figure  III-2  illustrated 


the  phenomena.  An  explanation,  first  proposed  by  Carl  Wunsch  (1977), 
which  is  consistent  with  this  latitudinal  structure  of  zonal  velocity 
as  well  as  with  the  steadiness  with  time  mentioned  as  point  1)  above, 
is  the  presence  of  a  first  meridional  mode  long  Rossby  wave. 

Only  two  and  a  half  weeks  elapsed  between  the  first  and  the  last  of 
these  six  transects,  so  that  it  is  not  surprising  that  no  significant 
phase  differences  at  depth  with  time  can  be  detected  by  examining  the 
latitudinal  sections.  Consequently,  the  periodicity  of  a  vertically 
standing  or  propagating  wave  cannot  be  determined  directly  from  the 
data.  However,  the  strongest  zonal  velocity  signal  in  the  ~400-day 
series  of  current  meter  data  from  the  western  equatorial  Indian  Ocean 
in  1979-1980  is  for  a  180-day  period  oscillation  (Luyten  and  Roemmich, 
1982).  That  period  will  be  used  as  an  estimate  to  examine  consistency 
relations,  since  the  1976  data  indicates  that  the  oscillation  is  of 
much  longer  period  than  the  length  of  this  time  series,  and  since  a  lot 
of  energy  is  input  at  that  frequency  into  the  Indian  Ocean  (Knox,  1976; 
McPhaden,  1982).  Note  that  if  a  periodic  function  is  responsible  for 
the  observed  signal,  then  increasing  the  period  length  beyond  the 
semi-annual  which  will  be  used  here  would  increase  the  zonal 
wavelength,  according  to  linear  theory,  but  would  have  little  impact 
upon  either  the  turning  latitude  or  the  radius  of  deformation,  and  thus 
would  have  little  effect  upon  the  latitudinal  structure  of  the  zonal 
velocity  field. 

Before  a  detailed  interpretation,  the  visual  and  statistical 
evidence  for  the  presence  of  strong  antisymmetry  will  be  discussed  more 
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thoroughly  than  in  Section  1 1 1— A.  All  six  of  the  latitudinal  transects 
extending  to  3*N  display  the  same  antisymmetry  between  the  profiles  at 
1.5*N  and  3eN.  With  vertical  profiles,  which  constitute  snapshots  of 
the  water  column,  one  is  limited  by  the  number  of  realizations  possible 
in  the  vertical.  Using  spectral  techniques,  for  example,  one  cannot 
discern  the  differences  in  energy  between  the  longest  vertical  scales 
(lowest  vertical  modes). 

The  scale  of  the  zonal  signal  in  this  data  is  a  significant  portion 
of  the  observed  water  column,  so  that  most  standard  statistical 
approaches  are  not  going  to  be  reliable  in  quantifying  the  picture.  As 
a  result,  the  measure  P  was  devised  in  Section  1 1 1— A,  and  used  to 
examine  the  strength  of  the  similarities  and  dissimilarities  between 
pairs  of  profiles,  by  summing  them  at  each  depth.  The  equatorial 
profiles  were  included  there  only  for  completeness,  since  the  vertical 
structure  at  the  equator  is  visibly  more  complex  than  that  at  3*N,  for 
example,  a  fact  «onsis.ent  with  the  general  expectation,  from 
equatorial  theory,  of  the  nature  of  equatorial  trapping,  and  also 
consistent  with  the  spectral  estimates  presented  by  latitude  in  Figure 
III-l. 

No  similar  antisymmetry  is  seen  in  the  series  of  Wiite  Horse 
profiles  collected  in  the  western  equatorial  Pacific  Ocean  in  1978 
(Eriksen,  1981).  The  latitudinal  distribution  and  level  of  zonal 
kinetic  energy  observed  in  the  two  series  are  indistinguishable, 
however,  and  a  fuller  comparison  of  the  two  sets  of  observations  is 
found  in  Chapter  V. 
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Table  III —4  sets  forth  some  of  the  relevant  consistency  checks 
possible  between  the  data  and  a  first  meridional  mode  long  Rossby  wave 
of,  for  tils  comparison,  six-month  period.  The  primary  points  of 
comparison  are  thus  the  zonal  scale  of  the  phenomenon,  the  latitudes 
where  a  180*  phase  reversal  occurs,  and  the  relative  amplitudes  of  the 
features  at  various  latitudinal  locations.  The  zonal  scale  cannot  be 
resolved  from  the  1976  data:  Figure  1 1 1—6  shows  pairs  of  near- 
simultaneous  profiles  from  the  two  equatorial  stations,  at  50°30'E  and 
53*E.  Two  weeks  elapsed  between  each  pair,  and  there  are  no  clear 
phase  differences  with  longitude,  or  with  time.  The  former  indicates 
that  the  zonal  scale  is  much  longer  than  the  separation  between  the 
stations  and  the  latter  is  consistent  with  a  period,  if  any,  of  at 
least  six  months,  i.e.,  six  times  the  elapsed  time  between  the  first 
and  last  of  the  profiles  at  each  station. 

The  second  poinc  of  comparison  is  the  latitude  of  phase  reversal, 
which  in  the  data  occurs  somewhere  between  the  stations  at  1.5*N  and 
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FIGURE  III —6 

Comparison  of  profiles  at  the  equator,  50*30' E  and  5: 


2*N. 

The  third  point  of  comparison  is  the  amplitude  decay  of  the  zonal 
velocity  associated  with  a  first  mode  Rossby  wave  from  its  primary 
maximum  at  the  equator.  Reference  to  Table  1 1 1—4  indicates  that  this 
is  where  the  major  discrepancy  in  this  interpretation  arises. 

Amplitudes  at  the  equator  and  3/4*  away  would  be  very  similar,  as  would 
amplitudes  at  1.5°  and  3*.  particularly  for  a  vertical  wavelength  of 
I500sm.  The  data  does  not  indicate  thi  jump  in  relative  amplitude 
between  these  pairs  of  latitudes,  however.  A  90  %  decrease  in 
amplitude  between  the  equator  and  3*N  would  incline  one  to  the  view 
that  a  feature  observable  at  the  equator  would  not  be  manifest  at  3*N. 

It  thus  stems  unlikely  that  the  antisymmetry  across  2* N  is  due  to 
the  presence  of  a  linear  first  meridional  mode  Rossby  wave.  The 
possibility  that  a  higher  odd  meridional  mode  Rossby  wave  is 
responsible  is  unlikely  because  a  node  in  the  zonal  velocity  field 
would  be  expected  between  3/4*N  and  1.5*N,  although  the  relative 
amplitudes  would  remain  higher  in  that  case.  Among  other  explanations 
are  the  possibility  that  the  WKB  approximation  is  becoming  less  valid 
for  the  longer  vertical  scales,  or  that  linear  theory  is  'tot  applicable 
to  this  particular  feature  of  the  vertical  structure.  In  addition,  the 
full  Impact  of  the  45*  slant  of  the  western  boundary  of  the  Indian 
Ocean  upon  the  velocity  field  is  not  known.  Moore  (personal 
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communication)  has  shown  that  a  tilt  of  the  western  boundary  alters  the 
symmetry  properties  of  the  reflection.  The  reflection  of  a  first  mode 
long  Rossby  wave,  for  example,  requires  a  mixed  Rossby-gravity  wave,  of 
opposite  meridional  symmetry,  to  satisfy  the  boundary  condition,  in 
addition  to  a  Kelvin  and  a  short  first  mode  Rossby  wave.  The  relative 
distribution  of  energy  shifts  in  favor  of  the  mixed  wave  as  the  tilt  of 
the  coast  to  the  east  increases. 
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CHAPTER  IV  -  Meridional  Velocity 


Introduction 

The  meridional  velocity  component  in  general  reveals  both  shorter 
vertical  length  scales  and  shorter  temporal  scales  of  variability  than 
the  zonal  component.  The  approach  in  this  chapter  will  be  to  describe 
and  discuss  information  obtained  from  the  White  Horse  profiles,  by 
means  of  standard  auto-  (Section  A)  and  cross-spectral  calculations 
(Section  B)  and  also  by  means  of  dropped  lagged  coherence  and  phase 
estimates  (Section  C).  Interpretation  of  these  results,  and  their 
evaluation  in  terms  of  linear  equatorial  wave  theory,  will  comprise 
Section  D. 

Section  A.  Latitudinal  Distribution  of  Meridional  Velocity 

A. 1,  Latitudinal  energy  density.  The  latitudinal  distribution  of 
meridional  kinetic  energy  density  over  various  vertical  wavenumber 
bands  is  condensed  into  Figure  I V— 1.  Each  point  represents 
autospectral  density  at  a  particular  latitude  averaged  over  five 
adjacent  wavenumber  bands.  A  concentration  of  energy  at  the  equator  is 
clearly  visible  at  all  but  the  lowest  wavenumber  estimates,  and  is 
particularly  strong  for  the  wavenumber  bands  centered  at  360  and 
240sm.  Note  that  there  were  only  three  profiles  made  at  3/4*S,  so  that 
the  spectral  estimates  for  that  latitude  have  larger  error  bars  and 
less  statistical  reliability  than  the  other  near-equatorial  locations. 
The  estimate  at  1. 5*N  is  statistically  lower  than  that  at  the  equator, 
at  the  95^  confidence  level,  for  the  720,  360,  and  240sm  bands,  and 


FIGURE  IV-1 


Autospectra  of  stretched  and  normalized  meridional 
velocity,  averaged  over  five  wavenumber  bands,  presented 
as  a  function  of  latitude,  from  3/4‘S  to  5*N. 
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spectral  energy  density  seems  to  level  off  to  a  stable  "background" 
level  at  that  latitude.  The  energy  concentration  and  its  scales  will 
be  examined  analytically  in  Section  0. 

The  behavior  of  the  180sm  band,  where  energy  concentration  is 
apparently  less  intense,  may  be  an  indication  that  the  estimates  are 
approaching  the  noise  level  of  the  stretched  and  scaled  measurements. 

It  was  mentioned  in  Chapter  II  that  the  maximum  stretched  interval 
between  data  points  was  63sm  in  the  depth  interval  used  for  the 
spectral  calculations,  so  that  the  shortest  resolvable  wavelength  would 
be  twice  that.  However,  since  for  this  dataset  only  two  transponders 
generally  were  available  for  the  triangulation  calculation  (see  Chapter 
1 1  —A*)  •  the  noise  level  in  v  may  be  higher  than  the  wave  energy  at 
wavelengths  somewhat  longer  than  130sm.  In  addition,  the  adaptive 
filter  described  in  Chapter  II-C,  which  is  used  to  smooth  the  stretched 
profiles,  smooths  the  lower  part  of  the  water  column  over  as  much  as 
150sm. 

A. 2.  Latitudinal  profile  sectiois.  The  latitudinal  structure  of 
meridional  velocity  profiles  along  53°E  is  shown  in  the  time  series  of 
meridional  sections,  extending  from  5*N  to  3/4‘S,  that  forms  Appendix 
A.  The  variability  of  vertical  structure  with  latitude  which  was 
evident  in  the  above  autospectra  (Figure  I V— 1)  is  clear  in  the  profiles. 

A  transect  from  the  middle  part  of  the  time  series,  section  ^4,  is 
snown  in  Figure  IV-2.  Both  straight  (v  in  cm/sec)  and  stretched  (v*  in 
scm/sec)  velocity  profiles  are  presented.  Individual  features  can 
sometimes  be  tracked  between  the  equatorial  and  near-equatorial 
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profiles,  but  usually  not  as  far  as  l^N,  although  this  visual 
approach  can  be  deceptive.  The  decrease  in  vertical  structure  away 
from  the  equator, i.  e.,  the  increase  in  dominant  vertical  scale,  is 
very  clear. 

A. 3.  Variance-preserving  Net  Autospectra.  Figure  I V— 3  presents  the 
autospectral  estimates  of  Figure  IV-1  in  greater  vertical  wavenumber 
detail.  All  profiles  at  a  particular  latitude  were  averaged,  but  there 
was  no  averaging  over  wavenumbers.  The  format  is  that  known  as 
variance-preserving,  since  the  area  under  the  curve  is  directly  related 
to  the  energy  in  the  wavenumber  band.  The  ordinate  is  energy  density 
times  vertical  wavenumber,  plotted  against  log  wavenumber  along  the 
abscissa. 

A  general  decrease  in  energy  away  from  the  equator  is  evident  at 
all  but  the  lowest  wavenumbers.  A  peak  is  apparent  at  1200sm,  which 
gradually  decreases  in  strength  away  from  the  equator  to  1.5°N,  where 
it  is  distinguishable  primarily  because  of  the  abrupt  decrease  in  level 
at  900sm.  A  secondary  peak  at  600sm  at  3/4'N  and  S  is  brc ader  at  the 
equator,  where  there  is  approximately  the  same  energy  density  at 
720sm.  Troughs  appear  at  900sm  and  at  514sm  at  the  three  equatorial 
and  near-equatorial  stations,  giving  the  peak  between  them  more 
significance  than  if  it  occurred  at  a  single  latitude.  A  broad  peak  at 
400-450sm  occurs  at  the  same  three  stations.  The  other  strong  feature 
is  a  peak  at  257sm,  which  occurs  only  at  the  equator. 

Although  the  peak  at  1200sm  is  somewhat  misleading  in  that  only  at 
the  equator  does  this  band  contain  more  energy  than  the  bands  on  either 
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FIGURE  IV— 3 

Autospectra  of  stretched  and  normalized  meridional 
velocity  (v*),  with  no  wavenumber  averaging,  presented 
for  each  net  as  a  function  of  vertical  wavenumber. 
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side,  at  all  four  latitudes  (3/4*S  to  1.5*N)  the  1200sm  band  is  more 
energetic  than  the  adjacent  higher  wavenumber  band,  centered  at  900sm. 
Examination  of  the  zonal  kinetic  energy  spectra  (Figure  II 1-3)  reveals 
a  trough  rather  than  a  peak  at  1200sm,  at  the  equator,  3/'  *N  and 
l.5*N.  At  each  of  these  latitudes,  the  meridional  component  is 
slightly  more  energetic  than  the  zonal,  but  in  no  case  is  the 
difference  statistically  significant. 

The  meridional  spectral  energy  at  450sm  is  approximately  twice  the 
zonal  energy  at  both  the  equator  and  3/4 *N.  For  both  velocity 
components,  there  is  slightly  more  energy  in  the  400sm  band  than  in  the 
450sm  band,  which  is  also  the  case  for  zonal  kinetic  energy  at  3/4*N. 
However,  meridional  energy  is  significantly  greater  at  450sm  than  at 
514sm  at  both  latitudes. 

Section  B.  Cross-spectra 

Coherence  between  stretched  variables  is  low  in  general. 

Significant  coherence  is  found  predominantly  at  lower  wavenumbers. 
Appendix  C  shows  all  significant  coherences  at  each  latitude  between 
the  three  variables  (u,  v,  J  ),  with  no  averaging  in  wavenumber  space 
(Table  C-l).  "Significant"  is  used  here  and  in  the  following  to  mean 
significantly  different  from  zero  true  coherence  at  the  95  %  confidence 
level.  Error  bars  on  the  phase  estimates  indicate  the  95 ^  confidence 
1 imits. 

Relationships  which  can  be  gleaned  by  examining  coherence  at  the 
same  location  in  unaver aged  bands  are  presented  in  Table  IV— 1.  "Leads" 
is  used  in  the  table  and  the  text  following  to  mean  leads  by  »/2  with 
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LATITUDE 

TABLE  IV-1 

Cross-spectra  of  all  three  variables  (u*,  v*,T*),  presented  as 
functions  of  latitude  and  vertical  wavenumber,  with  no  wavenumber 
averaging.  All  significant  coherence  pairs  are  presented,  with 
their  phase  estimates.  The  same  Information  is  contained  In  Table  1 
of  Appendix  C,  and  this  figure  is  repeated  from  Figure  1 1 1-4. 
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depth  (within  the  error  bars  at  the  95%  confidence  level).  "Out  of 
phase"  means  a  phase  estimate  of  -180*. 

When  the  cross-spectral  estimates  are  averaged  over  "two" 
wavenumber  bands,  only  5  %  of  the  coherences  between  meridional 
velocity  and  the  other  two  variables  are  significant  at  the  95% 
confidence  level.  Since  a  95  %  confidence  level  means  that  5%  of  the 
estimates  of  zero  true  coherence  will  appear  significantly  different 
from  zero,  this  reduction  in  coherence  implies  that  the  coherent 
processes  have  narrow  vertical  wavenumber  bandwidths.  Table  2  in 
Appendix  C  shows  all  significant  coherence  pairs  between  the  three 
variables  when  the  estimates  are  averaged  over  "two"  adjacent 
wavenumber  bands,  equivalent  to  weights  of  1/4,  1/2,  and  1/4  about  the 
center  band,  so  that  adjacent  estimates  are  not  independent. 

The  following  discussions  will  focus  primarily  upon  cross-spectral 
relations  in  the  wavenumber  bands  where  dropped  lagged  coherence  (to  be 
discussed  in  Section  C)  has  shown  a  strong  indication  of  periodicity. 
Bands  adjacent  to  these  primary  ones  are  scanned  to  ste  how  broad  the 
relationships  are  in  wavenumber  space. 

B.2.  1200sm  wavelength  cross-spectra.  Examining  coherence  in  the 
1200sm  band  between  orthogonal  velocity  components  from  profiles  at  the 
same  location,  one  sees  significant  coherence  at  the  equator  in  the 
1200sm  band  and  also  in  the  adjacent:  900sm  band.  The  phase  relation  is 
not  stable,  being  indistinguishable  from  180*  at  1200sm  and  from  -90* 
at  900  sm,  in  the  latter  case  Indicating  u  leads  v  with  increasing 
depth.  At  3/4*N,  the  velocity  components  are  coherent  over  a  broad 
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wavenumber  range,  from  wavenumber  bands  centered  )etween  600  and  1800sm. 
At  1200sm  the  phase  estimate  is  close  to  -180*  (-139**30*),  while  in  the 
neighboring  bands  the  phase  is  indistinguishable  from  zero. 

Assuming  independence  of  the  profiles  in  the  subset  for  which 
displacement  at  the  equator  is  available,  the  only  significant  coherence 
for  the  1200sm  band  is  between  velocity  components;  neither  velocity 
component  is  coherent  with  displacement.  The  phase  estimate  between  u 
and  v  for  the  subset  is  170**27°,  indistinguishable  from  that  for  the 
full  set  of  profiles.  Coherence  in  adjacent  bands  is  limited  to  vertical 
displacement-zonal  velocity  at  900sm,  although  all  three  variables  are 
coherent  at  720sm,  with  displacement  leading  both  u  and  v. 

At  3/4°N  for  the  1200sm  band,  again  only  the  velocity  components  are 
coherent.  The  phase  estimate  for  the  displacement  subset  is  -147**29*, 
compared  to  -139**30°  when  all  drops  are  used.  There  is  an  indication 
(not  statistically  significant)  that  displacement  leads  zonal  velocity  in 
th's  band,  as  it  does  in  the  bands  to  each  side,  where  coherence  is 
significant.  [At  1800$m,  .82(90**29°);  at  1200sm,  .55(61**180°);  at 
900sm,  . 71(l!2**43*),  compared  to  the  coherence  cutoff  of  .67].  There  is 
no  other  coherent  relationship  between  variables  at  1800sm,  and  at  900sm 
the  other  phase  relations  between  variables  are  not  distinguishable  from 
zero,  although  coherence  is  significant. 

B.3.  450sm  wavelength  cross-spectra.  By  contrast  with  the  1200sm 
wavelength  behavior,  the  cross-spectra  between  velocity  components  show 
no  coherence  in  the  450sm  band,  at  either  the  equator  or  3/4*N.  The 
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adjacent  longer  wavelength  band  is  also  not  coherent,  but  at  both 
latitudes  the  400sm  band  shows  significant  coherence,  with  a  phase 
estimate  of  -133*22*  at  the  equator  and  -88**35*  at  3/4*N  (both 
indicating  u  leads  v). 

None  of  the  cross-spectra  between  velocity  components  and 
displacement  shows  coherence  at  450sm  at  either  the  equator  or  3/4°N. 
The  same  holds  for  514sm,  the  adjacent  higher  wavelength  band.  At 
400sm,  however,  the  coherence  and  jhase  estimates  between  velocity 
components  for  each  subset  agree  with  the  above.  In  addition,  at  3/4*N 
displacement  at  a  400sm  wavelength  is  coherent  with  both  velocity 
components.  Phase  estimates  indicate  displacement  leads  u  (96**33°) 
and  that  v  and  displacement  are  out  of  phase.  Displacement  is  not 
coherent  with  either  velocity  component  at  the  equator  in  the  400sm 
band,  although  the  coherence  amplitude  between  u  and  displacement  is 
sufficiently  high  (.71  versus  cutoff  of  .73)  to  give  a  phase  estimate 
significantly  different  from  zero,  and  consistent  with  that  calculated 
at  3/4*N  (83**54°). 


Section  C.  Dropped  Lagged  Coherence  (DLC)  Results 


C.l.  General.  Several  dominant  features  emerge  from  the  DLC 
analysis.  In  contrast  to  a  mid-latitude  spectral  model  of  internal 
waves  at  low  wavenumber  (GM76:  Desaubies,  1976),  it  is  clear  that  these 
equatorial  dropped  lagged  coherence  estimates  are  a  function  of 


vertical  wavenumber.  In  general,  pairs  of  the  zonal  velocity  component 
show  high  coherence  for  the  longest  resolvable  vertical  wavelength 


bands  at  all  three  latitudes  for  which  there  are  enough  velocity 
profiles  to  use  the  temporal  lag  approach:  at  the  equator,  3/4*N  and 
1.5*N,  53*E.  For  the  meridional  component,  the  opposite  is  true. 
Coherence  is  reduced  for  the  lowest  vertical  wavenumbers,  and  in 
general  is  not  significant.  Consequently,  the  high  coherence  for  the 
meridional  component  in  the  1200  stretched  meter  band  is  prominent,  as 
is  shown  in  Figure  I V— 4.  Examination  of  the  corresponding  phase 
estimates  for  evidence  of  periodicity  reveals  a  strong  signal  in  the 
1200sm  band  for  the  meridional  velocity.  There  is  in  general  more 
evidence  of  periodicity  for  the  meridional  component  than  for  the 
zonal,  and  stronger  evidence  of  periodicity  at  the  equator  than  away 
from  it. 

These  general  results  are  consistent  with  first  impressions  from 
examination  of  the  original  unstretched  data:  a  shorter-term  temporal 
variability  evident  in  the  north-south  velocity  component  in  sharp 
contrast  to  the  persistence  of  features  in  the  east-west  component 
during  the  time  span  covered  by  the  profiles,  and  a  tendency  toward 
shorter  vertical  scales  in  the  meridional  velocity  component  compared 
to  the  zonal. 

C.2.  DLC  at  12Q0sm.  At  a  vertical  wavelength  of  1200sm,  the 
meridional  components  are  coherent  at  the  equator  and  at  3/4*N.  The 
adjacent  bands  are  not  coherent.  Steadily  decreasing  phase  with  longer 
temporal  lag  is  observed,  indicating  phase  propagation  downward.  (See 
the  discussion  of  Figure  1 1-9  for  the  interpretation  of  DLC  phase  in 
terms  of  propagating  waves.)  At  the  equator,  the  meridional  component 
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FIGURE  IV— 4 . 

Superposition  of  dropped  lagged  coherence  and#phase  for 
the  meridional  velocity  component  at  net  6,  0  ,  53  E. 

All  eight  resolvable  temporal  lags  are  presented,  from 
1  day  to  19  days.  The  consistent  high  coherence  for  the 
1200sm  band  Is  clear  for  the  first  six  lags;  the  steady 
decrease  in  the  phase  estimates  Is  apparent  through  the 
longest  lag.  The  darkest  curves  are  for  the  shortest 
lags,  and  vice  versa.  The  same  information  Is  presented 
In  different  form  for  the  1200sm  band  In  Figure  IV-5;  for 
450sm  In  Figure  IV-6;  and  for  720sm  In  Figure  IV-7. 


0°  45'  N,  53 


functions  of  temporal  lag,  for  the  1200sm  vertical 
wavelength  band.  The  cone  which  defines  the  inferred 
range  of  periodicity  is  indicated  by  cross-hatching. 
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is  coherent  at  six  of  the  eight  temporal  lags,  with  three  of  the  phase 
estimates  non-zero  at  the  95^  confidence  level.  At  3/4*N,  five  of  six 
lags  are  coherent,  with  three  of  the  phase  estimates  .ignif icantly 
different  from  zero.  See  Figure  I V— 5. 

Using  95  *  error  bars  on  meridional  phase  estimates  to  compute 
period  ranges,  at  the  equator  one  obtains  a  range  of  60—77  days;  at 
3/4*N,  57-161  days.  A  linear  least-squares  fit  to  the  phase  estimates 
at  both  latitudes  yields  a  joint  estimate  of  72  days,  whether  all  phase 
estimates  are  used  or  only  those  for  which  the  coherence  was 
significant. 

Dropped  lagged  coherence  was  also  computed  at  1.5’N,  where  only 
four  temporal  lags  were  resolvable.  The  meridional  velocity  component 
is  coherent  through  ten  days  there,  and  the  phase  estimates  are 
increasingly  negative,  although  not  significantly  difierent  from  zero. 

It  should  be  noted  that  at  all  three  latitudes  (0°,  3/4°N,  1.5°N) 
where  dropped  lagged  coherences  could  be  computed,  the  phase  estimate 
for  the  meridional  component  in  the  1200sm  band  contirues  to  decrease 
at  a  steady  rate,  even  when  the  coherence  estimate  is  no  longer 
significantly  different  from  zero. 

The  best  estimate  for  the  periodic  signal  in  the  DLC  results  for  a 
wavelength  of  1200sm  is  72  days,  from  the  above.  For  reference,  recall 
that  the  minimum  Rossby  wave  period  at  this  vertical  wavelength  is 
86  days;  the  maximum  inertia-gravity  period  is  15  days.  (These  limits 
were  derived  in  Chapter  I-B.2,  for  Table  1-2).  For  tha  second 
meridional  mode,  the  minimum  Rossby  wave  period  is  112  days,  maximum 
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inertia-gravity  11.  Only  Kelvin  and  mixed  Rossby-gravity  waves  satisfy 
the  linear  dispersion  relation  between  the  maximum  inertia-gravity  and 
the  minimum  Rossby  wave  periods.  A  Kelvin  mode  is  not  relevant  here 
since  the  signal  appears  in  the  meridional  velocity  component,  at  and 
near  the  equator. 

C.3.  DLC  for  450sm.  Strong  evidence  of  periodicity  is  also  found 
in  the  dropped  lagged  coherences  for  the  450sm  wavelength  band.  At  the 
equator,  coherence  is  significantly  different  from  zero  for  aggregated 
pairs  of  meridional  velocity  through  temporal  lags  of  eleven  days;  for 
five  of  these  six  estimates,  the  phase  is  also  significantly  different 
from  zero.  For  the  next  longer  lag,  12.5  days,  coherence  is  slightly 
below  the  cutoff  but  the  phase  estimate  is  still  significantly 
different  from  zero  and  consistent  with  those  preceding  it.  The 
inferred  propagation  direction  is  the  reverse  of  that  for  the  1200sm 
band:  in  this  case,  phase  is  steadily  increasing  with  temporal  lag, 
indicating  phase  propagation  upwards,  with  the  second  of  each  pair 
leading  the  first.  The  meridional  OLC  at  3/4*N  yields  phase  estimates 
that  are  consistent  through  temporal  lags  of  9  days  with  those  at  the 
equator,  but  only  two  of  the  six  coherence  estimates  are  significant 
and  their  significance  is  borderline  (.81,  .79  for  cutoff  .79).  See 
Figure  I V— 6. 

When  95  %  error  bars  on  the  meridional  phase  estimates  are  used  to 
compute  a  period  range  for  the  equatorial  data,  the  result  is  47-72 
days.  A  linear  fit  to  the  phase  estimates  for  which  coherence  is 
significant  yields  a  period  estimate  of  57  days.  (For  reference,  the 
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minimum  Rossby  wave  period  for  this  vertical  wavelength  is  141  days; 
the  maximum  inertial -gravity  is  24  days.) 

The  adjacent  longer  wavelength  band,  centered  at  514sm,  is  less 
coherent  meridionally  at  the  equator  and  phase  estimates  do  not 
progress  uniformly.  The  phase  estimates  are  included  within  the  error 
bars  of  the  estimates  at  450sm,  however,  and  three  of  eight  are 
significantly  different  from  zero.  The  adjacent  shorter  band,  centered 
at  400sm,  shows  coherence  over  the  same  temporal  lags  as  the  450sm 
band,  but  no  phase  estimates  are  non-zero.  There  is  no  coherence 
meridionally  at  3/4*N  in  the  514sm  band,  and  although  some  coherence  is 
found  in  the  400sm  band,  phases  are  indistinguishable  from  zero. 

C.4.  DLC  for  720sm  band.  A  third  wavenumber  band  where  the 


meridional  dropped  lagged  coherences  indicate  some  periodicity  is 
centered  on  720sm.  Coherence  is  significant  for  half  the  eight  lag 
groups  at  the  equator  and  two  of  the  corresponding  phases  are 
significantly  different  from  zero,  so  that  a  period  range  of  31-44  days 
can  be  computed  from  the  95%  error  bars  on  the  phase  estimates.  The 
phase  estimates  are  consistently  negative,  indicating  phase  propagation 
downward.  Fcur  of  the  six  coherence  estimates  are  significant  at 
3/4*N.  The  phase  estimates  are  uniformly  negative,  decreasing 
regularly  with  increasing  lag,  and  again  two  of  the  phase  estimates  are 
significantly  non-zero.  They  indicate  a  somewhat  longer  period  of  80 
days,  with  the  range  being  52-254  days,  so  that  the  two  sets  of  period 
estimates  do  not  overlap  (Figure  IV-7). 
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Section  0.  Interpretation  and  Examination  in  Terms  of  Linear  Theory 

D.l.  Equatorial  trapping.  As  shown  in  Chapter  I,  the  latitudinal 
structure  of  the  dominant  velocity  component  of  the  two  gravest 
equator i ally-trapped  meridional  modes  is  a  Gaussian  centered  at  the 
equator.  For  a  Kelvin  wave  (n=-l),  the  only  velocity  component  is  the 
zonal;  for  a  mixed  Rossby-gravity  wave  (n=0),  meridional  velocity  has 
the  same  structure  as  u  for  a  Kelvin  and  the  off-equatorial  maxima  of 
zonal  velocity  (which  has  a  node  at  the  equator)  are  of  lower  amplitude 
than  the  meridional  component. 

The  latitude  [(Nq/b  |m|)^^]  at  which  the  kinetic  energy  of 
the  Gaussian  velocity  component  falls  to  one-third  its  equatorial  value 
is  less  than  3/4*  for  vertical  wavelengths  less  than  5l4sm  (see  Table 
1-1,  equator  al  Rossby  radius  of  deformation).  Consequently,  the 
latitudinal  structure  of  Kelvin  and  mixed  Rossby-gravity  waves  with 
vertical  wavelengths  less  than  500sm,  if  they  are  present,  cannot  be 
well-resolved  with  this  data  set.  Cross-spectral  phase  relations,  in 
particular,  are  likely  to  be  dominated  by  noise  at  short  vertical 
wavelengths,  since  the  equator  is  the  location  for  a  node  in  one  of  the 
velocity  components. 

In  light  of  the  above,  the  equatorial  concentration  of  meridional 
kinetic  energy  evident  in  Figure  I V— 1  can  be  examined  to  see  if  the 
Rossby  radius  is  the  appropriate  scaling,  that  is,  whether  the 
concentration  of  energy  is  consistent  with  the  presence  of  mixed 
Rossby-gravity  waves.  Table  IV-2  shows  a  comparison  of  energy  density 
off  the  equator  to  the  equatorial  values. 
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TABLE  IV-2. 

Equatorial  trapping  of  meridional  velocity 


lr 

(km) 

(sm) 

Autospectral  estimates 

0°  3/4°N 

Energy  ratios 

est.3/4°N  low  3/4°N  /3/4V 
est.O6  high  0°  ev  LR  ' 

adj .3/4°N 
adj. 0s 

98 

720 

129(97-180)  118(87-155) 

.91 

.48  .48 

(.24) 

69 

360 

51(38-70)  27(20-38) 

.54 

.29  .24 

(.19) 

56 

240 

19(14.5-26)  11(8.2-15) 

.57 

.31  .11 

(.18) 

1.5°N 

est. 1 .5°N 

low  1 . 5°N  ,1.5°v2  adj.l .5°N 

est.O6 

high  0°  e'  LR  ' 

adj.O 

91 

720 

57 

.44 

.06 

6S 

360 

8 

.15 

.003 

56 

240 

4.2 

.12 

.0002 

(  ): 

95%  confidence  limits 

adj.:  estimate  at  1.5°N  subtracted  from  numerator  and  denominator  before 
division. 
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Column  7  is  the  proportion  of  equatorial  energy  that  would  be 
present  at  3/4*  off  the  equator  if  the  Rossby  radius  is  the  controlling 
latitudinal  scale  factor.  Column  5  is  the  quotient  of  the  spectral 
estimates  at  3/4*f  and  S  and  0*,  and  column  6  shows  the  minimum  ratio, 
produced  when  the  lower  bound  on  the  estimates  at  3/4*  and  the  upper 
bound  at  0*  are  used.  The  upper  and  lower  bounds  are  the  confidence 
limits  at  95^.  Comparing  columns  5  and  7,  one  can  see  that  the  actual 
energy  decrease  is  slower  than  the  Gaussian,  although  when  the  extremes 
at  95  *  confidence  are  used  (column  6),  energy  in  the  720sm  band 
decreases  to  3/4*N  according  to  the  Gaussian.  Energy  in  the  360sm  band 
falls  off  somewhat  more  slowly,  and  in  the  240sm  band  is  almost  three 
times  as  large  as  the  Gaussian  model. 

The  lower  part  of  the  table  concerns  the  estimates  at  1.5*N  and  it 
is  clear  that  no  amount  of  manipulation  will  extract  agreement  with  a 
Gaussian  decrease  in  energy.  However,  if  the  apparent  leveling  off  at 
1.5*N  is  taken  into  account  and  the  spectral  estimates  at  that  latitude 
are  considered  as  background  levels,  they  can  be  subtracted  from  the 
equatorial  and  near-equatorial  estimates  (see  column  8).  Energy  at 
360sm  falls  off  well  within  the  predicted  model,  and  at  240sm  is  closer 
to  the  prediction.  The  latitudinal  distribution  of  meridional  kinetic 
energy  in  the  720,  360  and  240sm  wavelength  bands  is  thus  consistent 
with  the  presence  of  mixed  Rossby-gravity  waves  at  those  vertical 
wavelengths. 

The  spectral  estimates  discussed  above  have  been  averaged  over  five 
adjacent  wavenumbers,  but  standard  coherence  and  phase  (Section  B)  and 
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dropped  lagged  coherence  and  phase  (Section  C)  indicated  significant 
behavior  in  individual  wavenumber  bands.  In  particular,  it  would  be 
interesting  to  be  able  to  explain  the  structure  of  the  peaks  in  the 
autospectra  shown  in  Section  A.  These  attributes  of  the  data  will  be 
synthesized  in  the  following  sections. 

D.2.  Interpretation  of  spectral  behavior  for  the  1200 ;m  wavelength 
oand.  Since  the  indication  of  temporal  periodicity  at  a  vertical 
wavelength  of  1200sm  is  solely  in  the  meridional  component,  a  tentative 
identification  of  a  mixed  Rossby-gravity  wave  can  be  made,  based 
primarily  upon  tne  period  estimate  from  the  equator,  which  is  more 
tightly  defined  than  that  at  3/4 *N  and  for  which  no  other  even 
meridional  modes  are  possible  with  real  zonal  wavenumber.  Choosing  a 
)lane  in  the  three-dimensional  (m,  k,  w)  dispersion  relation  (Equation 
[-10)  based  on  this  vertical  wavenumber,  and  cutting  that  plane  at  this 
frequency  range,  one  obtains  a  zonal  wavelength  of  300-400km  when  the 
meridional  mode  number  is  zero.  See  Figure  1-6.  Since  the  signal  is 
in  v  at  the  equator,  only  even  meridional  modes  [n= 0,2,4,...]  need  be 
considered.  For  n  >  0,  the  dispersion  relation  yields  a  complex  zonal 
wavenumber,  deemed  to  be  unrealizable  considering  the  distance  of  the 
profile  station  from  the  African  coast  (D.  Moore,  personal 
communication).  The  period  estimate  from  3/4*N  overlaps  that  at  the 
equator,  but,  being  broader,  also  admits  the  possibility  of  the  second 
meridional  Rossby  mode.  The  e-folding  distance  for  mixed 
Rossby-gravity  wave  meridional  kinetic  energy  at  this  vertical 
wavelength  is  126km,  i.e.,  energy  decay  to  65 %  at  3/4*,  18  %  at  1.5*. 
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The  turning  latitude  for  a  second  mode  Rossby  wave  is  280km. 

[y-j-  *  ( (2n+l)NQ/e  m  )^.]  The  fact  that  the  period  signal  occurs  in 
the  meridional  velocity  component  at  the  equator  rules  out  a  Kelvin  or 
first  Rossby  mode  being  responsible. 

The  analysis  of  the  White  Horse  profiles  can  be  presented  as  an 
information  matrix.  The  interpretation  problem  is  then  reduced  to  the 
finding  of  a  coherent  solution  consistent  with  the  various  elements.  The 
matrix  components  are  the  autospectra  of  the  three  variables,  £*3u,  [2]Vj 
[3]  5  ;  the  dropped  lagged  cross-spectra  of  [4]u*u>  [5]v*7  (the  number 
of  displacement  profiles  is  insufficient  to  compute  DLC  tor  3  );  and  the 
more  standard  cross-spectra  of  l-7^u*3,  and  C8]v*3  . 

In  the  context  of  linear  equatorial  wave  theory,  a  testable 
hypothesis  generated  from  [5],  the  DLC  for  v,  is  the  presence  of  a 
mixed  Rossby-gravity  wave  of  1200sm  vertical  wavelength,  with  downward 
phase  propagation  and  slightly  longer  than  two-month  period.  At  the 
equator,  zero  zonal  velocity  and  vertical  displacement  are  associated 
with  such  a  wave  (i.e.,  both  variables  have  a  node  at  the  equator),  and 
to  the  north  of  the  equator,  v  and  3  would  be  in  phase,  leading  u  by 
«/2  with  depth.  Examination  of  [2],  meridional  kinetic  energy  at  a 
wavelength  of  1200sm,  as  a  function  of  latitude  indicates  the  e-folding 
distance  from  the  equator  is  broader  than  the  predicted  Gaussian. 

There  is  a  falloff,  however,  as  well  as  evidence  that  the  apparently 
slow  rate  might  be  attributable  to  a  high  background  level. 

Examination  of  [4],  DLC  for  u,  reveals  that  no  periodicity  can  be 
directly  ascribed  to  the  zonal  velocity  component:  the  high  coherence 
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and  near  zero  phase  estimates  are  compatible  with  the  presence  of  a 
long-period  standing  wave,  with  a  period  greater  than  80  days,  or  of  a 
propagating  wave  of  longer  period.  External  evidence  (from  current 
neters  set  in  1976  and  1979)  indicates  a  general  concentration  of  zonal 
kinetic  energy  at  periods  greater  than  two  months  (Luyten,  1982),  and, 
in  particular,  a  highly  energetic  vertically  propagating  signal  of 
semi-annual  period,  with  a  vertical  wavelength  of  the  order  of  the 
ocean  depth  (Luyten  and  Roemmich,  1982).  Attenuation  with  depth  due  to 
weak  reflections  from  the  north-south  boundaries  might  yield  an  energy 
distribution  over  the  longest  resolvable  vertical  scales,  when  analyzed 
by  means  of  vertical  profiles,  consistent  with  the  high  DIC  coherence 
[4]  as  well  as  the  energy  level  [1]  in  those  bands  (Luyter ,  personal 
communication). 

The  absence  of  coherence  between  vertical  displacement  and  the 
velocity  components  ([7]  and  [8])  off  the  equator  is  consistent  with 
the  fact  that  the  calculation  of  j  is  a  noisy  one  and  also  with  the 
relatively  low  amplitudes  of  u  and  J  for  a  mixed  Ross by-gravity  wave 
of  this  period.  Although  the  ratio  u2/v2  increases  with  increasing 
latitude  for  a  mixed  Rossby-gravity  wave  (5  %  at  3/4*,  21  %  at  1.5*, 
for  a  72-day  period  at  this  vertical  wavelength),  the  ratio 
v2/v2( 0*)  falls  so  quickly  (to  18  %  by  1.5*)  that  any  zonal 
velocity  signal  would  be  lost  in  the  noise.  Displacement  for  a  mixed 
Rossby-gravity  wave  follows  the  same  pattern  as  does  the  zjnal 
component . 
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The  cross-spectra  of  velocity  components  [6]  seem  to  provide  the 
major  discrepancy  for  this  interpretation,  since  they  indicate  that  u 
and  v  are  coherent  and  180*  out  of  phase  at  both  the  equator  and  3/4*N, 
rather  than  the  expected  lack  of  coherence  at  the  equator  and  90*  phase 
difference  off  the  equator.  Additional  information  with  respect  to  the 
average  phase  difference  between  the  components  seems  to  support  the 
decoupling  interpretation  of  zonal  energy  above,  however.  When  the 
complete  series  of  profiles  are  considered  in  shorter  groups  and 
cross-spectra  are  computed  (shown  below  in  Table  IV-3),  the  phase 
difference  is  seen  to  be  a  function  of  time. 


TABLE  IV-3:  UV  Cross-spectra 

equator  #  3/4*N 

coherence  phase[*J  coherence  phase[*J 

all  profiles  .84( .59)  161*20  ,75(.59)  -139*30 

early  .84(.79)  -139*40  u  leads  . 77( . 79)  -115*56  u  leads 

middle  .87(.73)  170*27  .82(.67)  -147*29  u  leads 

late  .98( . 79)  156*12  v  leads  .97( .79)  -174*15  u  leads 

(  ):  95  %  cutoff  level  for  true  coherence  significantly  different  from 
zero. 

*  :  95  %  confidence  interval  on  phase  estimates. 


Consider  the  nth  harmonic  of  two  series,  up  =  acos(2nknz  =  <*n)  and 
vn  =  bcos(2irknz  +  en).  ©n(=  en  -  an)  is  then  the  average  phase 
lead  of  vn  over  un  with  depth  and  ep  =  f(t)  implies  that  en  and/or  <*n 
are  some  functions  of  time.  Let  Bp  =  -2u)t+Cj,  ap  =  -wt-*^;  then 

1  2  it 

®n  sr  -wt+Cy  Over  an  elapsed  time  of  at  =  "jg  say,  one-tenth  the 
period  of  the  longer-period  wave,  the  average  phase  difference  will 
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change  by  30*  and  in  the  observed  direction,  if  the  wave  frequencies 
are  related  as  above. 

The  choice  of  frequencies  in  the  above  analysis  is  deliberate,  with 
one  frequency  twice  that  of  the  other.  The  1976  White  Horse  profiles, 
analyzed  by  means  of  dropped  lagged  cross-spectra,  indicate  a  temporal 
periodicity  in  the  meridional  velocity  component  of  two  to  three 
months.  The  1979-80  current  meter  measurements  indicate  a 
vertically-propagating  wave  in  the  zonal  velocity  component  with  a 
period  of  six  months.  Although  the  observed  coherence  amplitude  and 
phase  estimates  for  a  vertical  wavelength  of  1200sm  cannot  be  explained 
by  the  presence  of  a  single  mixed  f ossby-gravity  wave,  an  infinite  sum 
of  wave  solutions  is  not  required  as  an  explanation.  Two  wave  forms 
will  suffice. 

The  apparent  inconsistency  in  an  interpretation  of  the  1200sm 
signal  in  energy  and  period  as  indicating  the  presence  of  a  mixed 
Rossby-gravity  wave  may  instead  be  jvidence  of  the  simultaneous 
presence  of  a  longer-period  wave,  with  energy  predominantly  in  the 
zonal  velocity  component.  Such  an  ad  hoc  approach  neglects  some  major 
dynamical  questions,  however,  among  them  linearity  and  energy  sources. 

The  validity  of  the  assumption  of  linearity  must  be  considered. 
Neglect  of  the  advective  terms  in  the  momentum  equation  requires  that 
the  phase  speed  be  large  compared  to  both  the  mean  flow  and  the 
perturbation  particle  speeds.  At  a  vertical  wavelength  of  1200sm  for  a 
mixed  Rossby-gravity  wave,  the  linear  dispersion  relation  yields  a 
zonal  wavelength  of  -404km  (phase  propagation  westward)  at  a  period  of 
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60  days,  or  a  zonal  wavelength  of  -316km  at  a  period  of  72  days.  The 
respective  zonal  phase  velocities  are  -7.8  and  -5.1  cm/sec,  which  are 
to  be  compared  to  rms  amplitudes  from  the  spectral  decomposition  of  the 
velocity  profiles  of  ~3.0  cm/sec  at  the  equator  for  this  wavenumber 
band.  It  is  clear  from  the  ratio  u/c  that  nonlinear  effects  may  not  be 
completely  negligible.  It  is  not  clear  what  their  impact  might  be  on  a 
mixed  Rossby-gravity  wave.  Evaluation  of  the  mean  flow,  in  order  to 
compare  its  magnitude  with  these  phase  speeds,  is  difficult  to  do  from 
the  profiles.  Although  the  individual  profile  means  are  all  of  order 
1-2  stretched  centimeters  per  second,  critical  levels  may  exist  at 
various  depths.  These  questions  will  be  topics  for  future  research, 
but  will  not  be  pursued  further  here. 

A  third  consideration  involves  the  penetration  of  energy  to  depth. 
For  a  linear  Kelvin  or  mixed  Rossby-gravity  wave,  the  ratio  of  zonal  to 
ver-.ical  group  velocities  is  NQ/wave  frequency,  exactly,  where  NQ  is  the 
buoyancy  frequency  used  to  stretch  and  normalize  the  original  velocity 
profiles  and  thus  represents  the  ambient  density  structure.  For  a 
mixed  Rossby-gravity  wave  of  60-day  period,  this  ratio  yields 
propagation  of  one  and  a  half  kilometers  zonally  for  one  stretched 
meter  vertically.  For  longer  periods,  the  angle  from  the  horizontal  at 
which  energy  would  propagate  is  even  smaller. 

The  correct  interpretation  of  the  dropped  lagged  coherence  results 
is  that  phase  propagation  is  downward  (profiles  earlier  in  time  lead 
subsequent  profiles),  implying  a  deep  rather  than  surface  energy 
source.  Source  possibilities  include  reflection  from  a  smooth  bottom 


or  bottom  topographic  features,  since  reflection  from  the  African  coast 


requires  preservation  of  frequency,  and  no  westward-propa jating  (i.e., 
incoming  to  the  African  coast)  linear  equatorial  waves  have  real  zonal 
wavenumbers  for  this  frequency  and  vertical  structure. 

It  should  be  noted  that  the  inferred  zonal  energy  propagation 
direction,  i.e.,  eastward,  has  been  based  exclusively  upon  the 
assumption  of  a  mixed  Rossby-gravity  wave  from  the  period  fit.  There 
is  no  direct  evidence  regarding  zonal  structure  or  propagation. 

However,  the  zonal  antisymmetry  which  is  apparent  in  Figure  I V— 8  for  v 
and  v*  offers  indirect  evidence  of  a  zonal  wavelength  approximately 
twice  the  station  separation.  The  figure  shows  three  pairs  of  velocity 
profiles:  all  of  the  profiles  are  from  the  equator,  and  the  pairs  were 
collected  at  approximately  two  week  intervals.  The  first  profile  in 
each  pair  is  from  O',  50'30'E;  the  second  is  from  O’,  53*E,  so  that  the 
profiles  are  taken  from  longitudes  -275km  apart.  Statistical 
resolution  of  the  associated  vertical  length  scale  is  complicated  by 
the  fact  that  too  few  profiles  exist  from  the  station  nearest  the 
African  coast  for  sufficient  degrees  of  freedom.  A  vertical  length 
scale  of  much  less  than  1200sm  is  apparent,  which  may  be  superimposed 
upon  a  longer  vertical  scale  oscillation. 

D.3.  Interpretation  of  spectral  behavior  in  the  450sm  wavelength 
band.  The  DLC  results  also  showed  a  strong  signal  at  450sm.  In  the 
framework  of  linear  equatorial  wave  theory,  at  a  scale  of  450sm  only  a 
mixed  Rossby-gravity  wave  displays  the  observed  periodicity  with  real 
zonal  wavenumber,  and  signal  in  the  meridional  component.  Given  a 
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FIGURE  IV— 8 

Comparison  of  pairs  of  profiles  at  the  equator:  500 
series  are  from  50  30'E;  600  series  are  from  53*E.  Each 
pair  was  separated  by  elapsed  time  of  two  weeks  from  the 
adjacent  palr(s). 
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vertical  wavelength  of  450sm  and  vertical  phase  propagation  upward,  a 
mixed  wave  with  47-day  period  will  have  a  zonal  wavelength  of  almost 
'000km  and  an  eastward  zonal  group  velocity  of  5.9cm/sec  or 
S.llkm/day.  Thus  almost  four  wave  periods  would  be  required  to  travel 
from  the  African  coast  to  53°E.  At  that  longitude,  it  would  have 
penetrated  ~800sm  deeper  into  the  water  column  while  traveling  a 
thousand  kilometers  zonal ly. 

If  the  longer  period  estimate,  72  days,  is  taken  from  the  dropped 
lagged  coherence  results,  the  zonal  wavelength  is  much  shorter,  ~410km, 
the  group  velocities  are  almost  halved,  and  depth  penetration  at  53°E 
is  about  535sm. 

Choosing  the  intermediate  period  of  57  days  from  the  linear  fit 
yields  a  zonal  wavelength  of  ~730km  and  zonal  group  velocity  of 
5cm/sec.  A  disturbance  originating  at  the  African  coast  with  energy 
propagation  downward,  related  to  the  Somali  current,  for  example,  and 
with  an  initial  depth  penetration  of  400m,  say,  will  have  propagated 
approximately  1090km  by  the  time  it  reaches  53°E.  The  depth  interval 
from  the  surface  to  400m  is  equivalent  to  the  interval  from  the  surface 
to  I700sm  in  the  stretched  coordinate,  and  in  1000km  such  a  disturbance 
would  have  penetrated  only  700sm  deeper  into  the  ocean.  Since  the 
portion  of  the  stretched  ocean  that  has  been  examined  in  the  above 
spectral  analysis  is  1400-5000sm,  a  coastal  disturbance  generating  a 
Kelvin  or  mixed  wave  at  this  period  (57  days)  will  only  have  involved 
the  upper  lOOOsm  of  the  examined  water  column.  For  a  longer  period, 
the  angle  from  the  horizontal  of  the  ray  path  along  which  energy 
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travels  is  even  shallower.  Reexamination  of  Figure  IV-8,  shown  earlier 
in  che  discussion  of  the  1200sm  wavelength  band,  demonstrates  that  the 
clearest  antisymmetry  is  in  the  upper  part  of  the  water  column  and  at  a 
lenqth  scale  shorter  than  1200sm,  so  that  the  longitudinal  reversal  of 
phase  may  relate  to  this  450sm  oscillation. 

For  the  above  periods,  the  relationship  between  zonal  phase  speed 
and  rms  zonal  velocity  is  consistent  with  the  linearization  assumption: 
the  rms  zonal  velocity  estimated  from  the  autospectra  is  ~.3cm/sec  at 
the  equator  aid  3/4*N,  at  least  one  order  of  magnitude  smaller  than  the 
phase  velocities  (~7cm/sec  for  a  period  of  72  days,  ~45cm/sec  for  a 
period  of  47  days).  The  mean  flow  .is  also  of  much  smaller  magnitude 
than  these  phase  speeds. 

The  energy  e-folding  distance  [NQ/fi  jmi]1^2,  i.e.,  the  equatorial 
Rossby  radius  of  deformation,  is  77  km,  so  it  is  not  surprising  to  find 
that  the  period  signal  is  in  fact  weaker  at  3/4*N. 

Conclusion 

There  is  strong  evidence  of  vertical  propagation  from  the  dropped 
lagged  cross-spectra  of  1976  White  Horse  profiles  of  the  meridional 
velocity  component.  The  indicated  propagation  is  at  two  different 
vertical  wavelengths  and  is  in  opposite  directions.  The  above 
interpretations  showed  that  all  of  the  information  available  from  the 
data  studied  here  are  at  least  reasonably  consistent  with  the  presence 
of  two  separate  mixed  Rossby-gravity  waves. 
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Chapter  V.  Conclusion 

The  major  results  from  this  research  are  the  evidence  of 
equatorial  trapping,  which  varies  with  vertical  wavenumber,  and 
evidence  of  vertical  propagation,  which  in  some  cases  is  consistent 
with  the  presence  of  particular  equator! ally-trapped  meridional  wave 
modes. 

The  presence  of  deep  energy,  well  below  the  thermocl ine,  was 
revealed  by  the  White  Horse  profiles  in  1976.  The  velocity  structure 
which  became  visible  was  not  any  simple  current  field  but  a  complex 
multi-directional,  rather  than  zonal ly-dominated  field,  albeit  with 
multiple  reversals.  A  unique  forcing  regime  exists  at  the  equator, 
.uch  that  waves  can  be  forced  at  any  frequency,  so  that  coupling  with 
itmospheric  variability  at  all  periods  is  possible.  With  respect  to 
fhe  Indian  Ocean  as  a  dumping  ground  for  atmospheric  enercy,  Luyten 
and  Roemmich  (1982)  have  estimated  that  the  downward  flux  of  energy 
through  the  thermocl ine  at  180  days  corresponds  to  an  amoint  15% 
of  the  annual  mean  energy  put  into  the  North  Atlantic. 

The  questions  of  what  happens  to  that  energy  put  in  al  the 
surface,  presumably  by  wind  forcing,  and  of  the  setup  prot lem  of 
modes,  particularly  meridional  modes  in  the  equatorial  waveguide, 
have  not  been  addressed  here.  The  setup  of  vertical  modes  versus 
vertical  propagation  has  been  examined  in  passing  and  remains  to  be 
comprehensively  addressed.  The  problems  involved  here  are  clear: 


only  in  special  cases  can  energy  in  a  particular  meridional  mode 
propagate  even  predominantly  vertically,  so  that  zonal  diffusion  of 
energy  at  a  particular  vertical  wavelength  and  particular  period  is 
involved.  The  x,z-angle  at  which  energy  propagates  requires  multiple 
reflections  from  sidewalls  (of  non-longitudinal  orientation  in  the 
western  Indian  Ocean)  as  well  as  from  the  bottom,  before  energy  could 
return  to  the  same  x,z-location.  Related  issues  are  the  mechanisms  for 
vertical  scale  selection:  the  problems  of  energy  propagation  through 
the  thermocline  and  through  the  zonal  surface  and  near-surface  currents 
have  also  not  been  addressed  here. 

What  has  been  attempted  is  a  definition  of  what  was  going  on 
beneath  the  thermocline  in  the  western  Indian  Ocean  during  May  and  June 
of  1976,  with  the  presumption  that  that  period  was  not  atypical. 

Energy  was  clearly  present:  what  was  its  form?  Initially,  an  attempt 
to  examine  possible  forcing  mechanisms  was  planned.  The  equatorial 
concentration  of  energy  seems  to  rule  out  off-equatorial  sources, 
although  the  latitudinal  coverage  is  limited  to  3/4°  south  of  the 
equator.  Ray  path  tracing  from  sources  such  as  the  Somali  Current  is 
constrained  by  the  direction  of  propagation  (phase  downward)  for  the 
strongest  and  most  consistent  signal  of  vertical  propagation. 

The  clearest  evidence  for  vertical  propagation  is  from  the 
meridional  velocity  component  at  vertical  wavelengths  of  1200sm  and 
450sm.  In  both  cases,  the  indicated  periods  allow  only  mixed  Rossby- 
gravity  waves  in  a  linear  wave  model.  The  autospectral  peak  at  1200sm 
for  v  is  echoed  by  the  peak  in  dropped  lagged  coherence  at  that 
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wavelength.  The  period  indicated  is  2-3  months,  with  a  best  fit  of  72 
days.  Energy  decreases  away  from  the  equator,  at  a  slightly  lower  rate 
than  the  theoretical  Gaussian  with  latitude.  The  zonal  wavelength  of 
300-400km  obtained  via  the  dispersion  relation  is  reinforced  by  a 
comparison  of  profiles  separated  by  2.5°  of  longitude  at  the  equator. 
Phase  propagation  is  downward,  implying  upward  energy  propagation. 

The  major  inconsistency  in  the  1200sm  interpretation,  high 
coherence  between  u  and  v  at  the  equator  with  a  phase  estimate  of 
approximately  18U*,  can  be  explained  by  the  simultaneous  presence  of  a 
Kelvin  or  Rossby  wave  of  at  least  twice  the  period  of  the  mixed  wave, 
with  approximately  the  same  vertical  scale.  This  explanation  is 
consistent  with  analysis  of  a  slightly  reduced  portion  of  the  water 
column,  which  resolves  different  vertical  wavenumber  bands,  and 
indicates  that  the  zonal  energy  is  concentrated  closer  to  a  vertical 
wavelength  of  1600sm  than  1200sm  and  is  consequently  not  well-resolved 
by  the  discretization  from  spectral  analysis  of  the  entire 
sub-thermocline  water  column. 

The  second  strongest  evidence  of  vertical  propagation  is  for  a 
vertical  wavelength  centered  at  45(<sm,  and  the  indicated  direction  of 
propagation  is  upward  in  this  case,  implying  a  surface  energy  source. 
The  energy  e-folding  distance  is  consistent  with  the  theoretical 
structure  of  a  mixed  Rossby-gravity  wave  and  the  period  estimate  from 
the  dropped  lagged  coherence  results  is  1  1/2-2  1/2  months,  with  a  best 
fit  of  57  days.  The  corresponding  zonal  wavelengths  calculated  using 
the.  linear  dispersion  relation  are  410-2000km,  with  57  days  yielding  a 


zonal  wavelength  of  ~730km. 

Dropped  lagged  coherence  analysis  of  the  zonal  velocity  component 
suggests  processes  of  longer  period  than  can  be  resolved  with  this  set 
of  observations.  The  phase  reversal  across  2*N  suggests  the  presence 
of  a  first  meridional  mode  Rossby  wave.  However,  the  energy  observed 
at  3*N  relative  to  the  equator  is  higher  than  explainable  by  the 
presence  of  a  first  mode  Rossby  wave. 

Examination  of  spectral  energy  density  variation  with  latitude 
showed  equatorial  intensification  at  all  vertical  wavenumbers  for  both 
velocity  components.  The  meridional  velocity  component  settles  to  an 
apparent  background  level  at  1.5*N,  and  a  Gaussian  decrease  in  energy 
away  from  the  equator  can  be  extracted  for  the  wavenumber  bands 
encompassing  287sm  through  1440sm.  That  decrease  is  consistent  with 
the  dominance  of  the  meridional  velocity  field  by  mixed  Rossby-gravity 
waves.  Note  that  the  slant  of  the  east  African  coast  near  the  equator, 
approximately  45*  to  the  meridian,  theoretically  favors  the  selection 
of  mixed  Rossby-gravity  waves  when  incoming  energy  is  reflected  (Moore, 
personal  communication). 

Zonal  kinetic  energy  is  less  tightly  trapped  to  the  equator  than 
meridional,  and  appears  to  decrease  more  gradually  to  5*N,  the  further¬ 
most  data  station.  For  all  four  of  the  averaged  wavenumber  bands 
studied  here,  however,  the  estimate  at  1.5*N  is  significantly  lower 
than  that  at  the  equator,  at  the  95%  confidence  level.  There  is 
relatively  more  energy  in  the  zonal  component  at  longer  wavelengths, 
relatively  more  in  the  meridional  component  at  shorter  wavelengths. 
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The  equatorial  intensification  is  stronger  at  shorter  wavelengths  for 
both  components,  consistent  with  linear  theory.  A  general  conclusion 
from  the  analysis  of  equatorial  trapping  is  that  more  than  just  mixed 
Rossby-gravity  and  Kelvin  waves  are  present.  The  distribution  of 
energy  implies  the  presence  of  Rossby  waves  as  well. 

A  comparable  series  of  profiles,  also  using  the  White  Horse,  was 
collected  in  the  western  equatorial  Pacific  Ocean  in  1978  by  Charles 
Eriksen  (1981).  That  series  has  greater  longitudinal  coverage,  so  that 
individual  stations  were  not  repeated  as  often  as  in  the  Indian  Ocean. 
Consequently,  those  profiles  cannot  be  examined  as  a  time  series.  The 
Pacific  profiles  cover  approximately  the  same  latitudinal  range  about 
the  equator,  so  that  comparisons  can  be  made  on  that  basis. 

Equatorial  intensif ication  was  observed  in  the  Pacific  data  for  all 
three  variables,  both  of  the  horizontal  velocity  components  and 
vertical  displacement.  The  relative  intensity  of  trapping  was  the 
reverse  oi  that  observed  in  the  Indian  Ocean,  i.e.,  stronger  for  the 
zonal  velocity  component  than  for  the  meridional.  The  displacement 
spectra  in  the  Indian  Ocean  do  not  exhibit  equatorial  intensification. 
The  averaged  autospectra  of  zonal  kinetic  energy  from  the  Pacific  and 
the  Indian  Oceans  are  not  distinguishable  either  in  energy  level  or  in 
latitudinal  distribution,  within  the  95^  confidence  level  error  bars 
(cf.  Eriksen,  1981,  figure  7,  and  Figure  III-l).  The  Indian  Ocean 
autospectral  estimates  for  meridional  kinetic  energy,  however,  are 
significantly  higher  at  the  equator  for  the  720,  360,  and  240sm 
wavelength  bands,  although  indistinguishable  from  the  Pacific 
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observations  1.5*  away  from  the  equator. 

A  major  result  of  Eriksen's  analysis  (1981,  1982)  was  the 
conclusion  from  the  phase  relation  between  zonal  velocity  and  vertical 
displacement  at  the  equator  (where  ^  led  u  by  90*  with  depth  for 
vertical  wavelengths  longer  than  lOOsdb)  and  from  the  equatorial 
trapping  of  those  two  variables  that  Kelvin  waves  dominated  near  the 
equator.  No  such  conclusion  can  be  drawn  for  the  Indian  Ocean  from 
this  data. 

Heavy  averaging  over  vertical  wavenumbers  was  required  for  the 
Pacific  data  in  order  to  obtain  statistically  significant  results, 
since  there  were  fewer  profiles  at  a  particular  station.  The  resultant 
spectra  are  smooth  and  significant  coherence  estimates  are  rare.  The 
discrete  vertical  structures  revealed  in  the  Indian  Ocean  spectra 
disappear  with  wavenumber  averaging,  as  does  the  significance  of  the 
coherence  pairs,  so  that  it  is  clear  that  the  dominant  vertical 
processes  in  the  Indian  Ocean  are  narrow  band  in  vertical  wavenumber. 


Our  observational  capacity  in  the  ocean  is  limited.  Current  meters 
can  be  deployed  for  long  time  spans  and  provide  excellent  sources  for 
monitoring  periodic  oscillations  at  specific  depths,  but  resources 
limit  those  specific  depths.  Vertical  profilers  offer  snapshots  of  the 
water  column,  but  using  them  to  determine  periodicity  or  to  examine  the 
longest  vertical  scales,  corresponding  for  example  to  the  lowest 
vertical  modes,  is  extremely  difficult.  Dropped  lagged  coherence,  used 
here,  seems  the  only  option  for  the  former  and  the  latter  is 
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complicated  by  the  inapplicability  or  gradual  breakdown  of  the  WKB 
approximation  as  the  scales  approach  those  of  the  vertical  variability 
of  the  buoyancy  frequency. 

In  any  case,  profilers  are  limited  by  the  number  of  realizations 
possible  between  the  surface  and  the  bottom  of  the  ocean.  Current 
meters  merely  need  to  be  deployed  some  multiple  of  the  period  of 
interest  in  order  to  obtain  multiple  realizations  in  time.  On  the 
other  hand,  one  of  the  problems  revealed  by  these  particular  profiles 
is  that  motion  at  depth  appears  highly  organized,  whereas  linear  theory 
indicates  that  energy  propagation  occurs  at  shallow  angles  to  the 
horizontal  for  low  meridional  modes.  Intuitively,  it  is  easier  to 
envision  the  setup  of  the  lower  modes. 

Additional  observations  exist  now  of  the  complex  vertical  structure 
in  the  equatorial  oceans.  The  particular  instrument  used  here  has  also 
been  used  in  the  equatorial  Pacific  Ocean  (Eriksen,  1981),  and  other 
profilers  have  been  used  in  both  the  Pacific  and  the  Indian  Oceans 
(Hayes,  1981;  Leetmaa  and  Spain,  1981).  Other  investigations  have 
discovered  the  dichotomy  in  frequency  dominance  of  the  velocity 
components.  Weisburg  and  Horigan  (1981)  discuss  evidence  in  the 
eastern  equatorial  Atlantic  for  a  vertically-propagating  mixed  Rossby- 
gravity  wave  of  30 -day  period  dominating  the  meridional  velocity 
component,  while  the  zonal  component  indicates  a  periodicity  on  the 
order  of  four  to  eight  months.  With  other  types  of  instruments, 
oscillations  with  periodicities  longer  than  several  weeks  have  been 
discovered  in  the  Indian  Ocean;  a  30-60  day  period  oscillation,  perhaps 
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identifiable  as  a  Kelvin  wave,  and  a  periodicity  of  180  days  (McPhaden, 
1982),  the  last  of  which  agrees  well  with  the  intense  zonal  signal 
analyzed  as  a  combination  of  Kelvin  and  first  mode  Rossby  wave  by 
Luyten  and  Roemmich  (1982).  The  latter  also  see  a  strong  signal  in  the 
meridional  velocity  component,  identified  as  a  mixed  Rossby-gravity 
wave  of  26-day  period  (Roemmich  and  Luyten,  personal  communication). 

Evidence  of  vertical  wave  propagation  in  the  deep  equatorial  oceans 
leads  to  questions  that  go  beyond  linear  equatorial  wave  theory.  The 
basic  state  of  the  equatorial  oceans  is  not  quiescent,  but  includes 
both  vertical  and  meridional  shear  flows.  Their  consequences  have  been 
well  studied  in  the  atmosphere,  but  the  oceanic  case  has  received  much 
less  attention.  Since  the  equations  of  motion  are  not  separable  for 
the  case  of  mean  vertical  shear,  a  numerical  model  is  required. 

The  analyses  of  wave  propagation  in  a  mean  flow  field  that  have 
been  done  to  date  for  the  equatorial  oceans  have  been  primarily 
concerned  with  latitudinal  variations  of  the  mean  flow  in  a  two-layer 
ocean  (McPhaden  and  Knox,  1979;  Philander,  1979).  Some  qualitative 
analysis  has  been  done  by  Philander  (1979)  of  the  effects  of  vertical 
shear  across  an  interface,  but  specific  examples  need  to  be  analyzed. 
Philander  (1978)  has  also  examined  the  transmission  of  equatorial  .waves 
through  a  typical  tropical  density  profile,  and  concludes  that  only 
particular  wavelengths  can  penetrate  to  great  depths.  Spectral 
analysis  indicates  certain  preferred  scales  in  the  western  Indian 
Ocean,  so  that  the  selection  process  becomes  a  first-order  question. 
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APPENDICES 


Appendix  A:  Eight  longitudinal  transects  showing  zonal  and 

meridional  velocity  profiles,  vertical  displacement 
profiles  in  original  form  and  after  WKB  stretching  and 
normalizing. 

1.  u,  v,  ^  from  5°N  to  3/4°S 

2.  u*.  v*,  -J*  from  5°N  to  3/4°S 


Appendix  B:  All  profiles  at  each  profiling  site  (i.e,  temporal 

transects  for  each  latitude): 

1.  u,  v,  J  for  each  net 

2.  u*.  v*,  for  each  net 


Appendix  C:  Tables  of  coherence  and  phase 

Table  1.  (no  wavenumber  averaging) 

a.  u*  *  v*  [all  profiles] 

b.  *  u*  [subset  for  which 

c.  *  v*  displacement  is  available] 

Table  2.  (averaged  over  2  adjacent  wavenumber  bands) 
u*  *  v*,  “J*  *  u*,  "J*  *  v*  [subset] 
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Eight  latitudinal  transects  showing  zonal  and  meridional  velocity 
and  vertical  displacement  profiles  in  both  original  form  and 
after  WKB  stretching  and  normalizing. 

For  each  section, 

a)  u,  v,  from  5*N  to  3/4°S. 

b)  u*,  v*,  from  5*N  to  3/4*S. 

The  profiles  at  the  end  points  are  common  tc  sequential 
sections.  The  distribution  of  the  sections  in  time  is  shown 
in  the  first  figure  (cruise  track).  Maximum  elapsed  time 
between  profiles  presented  in  the  same  section  was  3.6  days 
(3/4°S  to  5”N),  and  average  elapsed  time  between  the  most 
closely  spaced  profiles  (3/4*S  to  1.5*N,  latitudinal  spacing 
3/4°)  was  12  hours. 
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APPENDIX  B 


All  profiles  at  each  profiling  site  (i.e.,  temporal  transects  for 
each  latitude): 

a)  u,  v  }  $  for  each  net 

b)  u*,  v*,  J*  for  each  net 
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APPENDIX  C 


Tables  of  Coherence  Pairs 
At  Each  Latitude 


Table 

a) 

b) 

c) 


1.  No  wavenumber  averaging 
u*v  fall  profiles] 

[*u  [subset  for  which  vertical 

;*v  displacement  is  available.] 


U’ 

r 


Table  2.  Averaged  over  two  adjacent  wavenumber  bands 
a)  u*v 


b) 

c) 


1*u 

1*' 
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TABLE  C-la 


u*v 


3/4*S 

0* 

3/4*N  1 .5*N 

3*N 

5°N 

Coh. (Phase) 

3600 

1800 

1200 

900 

720 

. 75( —40 ) 

.84(166*22) 

.63(57*45) 

.65(-34±42) 

.75(-134±30) 

.69(-28±37)  . 76(— 2  3 ) 

.75(34*30) 

.  76( — 11 ) 

.74(37) 

.88(5) 

.89(— 6 ) 
.94(17) 
.84(-34) 

600 

514 

450 

400 

360 

None 

.83( -133*23) 

.74(-54*31) 

.71 (-88*35) 

327 

300 

277 

257 

240 

.60(119*50) 

.58(164*54) 

.60 (-34*50) 

.76(110) 

.74(109) 

225 

212 

200 

189 

180 

171 

164 

. 65( —13*42 ) 

OOOOO  0  0-000  r-^  O  r--  r-~  O  in  M  O  O'  O  o  ^ 

OOOON  ooino®  n  O  ivirnf  (v  o  o  03  oo  rv  in 

co  co  cnj  O' iOui<tOf o  ro  ro  cm  cm  cm  cm  cm  cm  .— t .— 1 .— t 


TABLE  C-lb 


}*u 


3/4*S  0’  3/4*N  1 .5*N 


*z  Coh. (Phase) 


.82(-9Q*31) 

.90(-148*23)  -72( —122*44) 
.89(-l29*25)  - 92(— 173)  .79(134*33) 

•69(-56*48) 


1.00(168)  .79( -96*33) 

. 98( —152)  ,77(-120*36) 


,83(-114*28) 

.99(159)  .82( -131*35)  .77(-56*36) 


.87 (-89*24) 

.85( -143*30) 

.84(-130*32)  .71(167) 


3*N 

.74(8) 

.94( -43*17) 
.74(178*48) 

.76(151*44) 

,76(-81*44) 


.73(50) 
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TABLE  C-lc 


1*V 


3/4*S 

£ 

3/4  *  N  1 . 5  *N 

3*N  5*N 

Coh. (Phase) 

3600 

.98(109) 

1800 

.82(-9*35) 

1200 

900 

•70(-146*46) 

720  .80 (-71*38)  .81(-130±31) 


600 

514 

450 

400 

360 


•84(-161*27) 


.98(123) 


327 

300  .76(46*36) 

277 

257 

240 


225 

212  .76(-32*44) 

200  .70 (-2 2*46) 

189  .70(28) 

180  .87( —5*27 ) 

171  .8l(-102*36)  . 94 ( — 136 ) 
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TABLE  C-2.  Coherence  Pairs  At  Each  Latitude 
Averaged  Over  Two  Wavenumber  Bands. 


3/4*S 

Coh. (Phase) 


3/4°N 
]  *u 


1.5*N 


1800 

.52  (-52) 

.72 

900 

.84( —142 ) 

•61( -141) 

600 

.63( -169) 

.63 

450 

.50 (-106) 

360 

•85(-173) 

. 71 ( —114) 

300 

.81(169) 

257 

.57 ( — 139 ) 

225 

. 5  4 ( —9 1 ) 

200 

.  59( — 163 ) 

180 

164 

.67 ( —140) 

.50(166) 

T*v 

1800 

J 

900 

. 57(— 150) 

.69 

600 

450 

.53(146) 

360  None 

300 

257 

225 

300 

180 

164 


•50(-105)  .76 
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